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ABSTRACT 

We obtain cosmological constraints from a measurement of the spherically aver- 
aged power spectrum of the distribution of about 90000 luminous red galaxies (LRGs) 
across 7646 deg^ in the Northern Galactic Cap from the seventh data release of the 
Sloan Digital Sky Survey. The errors and mode correlations are estimated thanks to 
the 160 LasDamas mock catalogues, created in order to simulate the same galaxies 
and to have the same selection as the data. We apply a model, that can accurately 
describe the full shape of the power spectrum with the use of a small number of free 
parameters. Using the LRG power spectrum, in combination with the latest mea- 
surement of the temperature and polarisation anisotropy in the cosmic microwave 
background (CMB), the luminosity-distance relation from the largest available type 
la supernovae (SNIa) dataset and a precise determination of the local Hubble pa- 
rameter, we obtain cosmological constraints for five different parameter spaces. When 
all the four experiments are combined, the flat ACDM model is characterised by 
^M = 0.259+°:ol5, ^b = 0.045 ± 0.001, n^ = 0.963 ± 0.011, ag = 0.802 ± 0.021 and 
Hq = 71.2 ± 1.4kms~^ Mpc~^. When we consider curvature as a free parameter, we 
do not detect deviations from flatness: 51k = (1-6 ± 5.4) x 10^^, when only CMB and 
the LRG power spectrum are used; the inclusion of the other two experiments do not 
improve substantially this result. We also test for possible deviations from the cosmo- 
logical constant paradigm. Considering the dark energy equation of state parameter 
Wde as time independent, we measure wue = ~l-025lQgg5, if the geometry is assumed 

to be flat, wde = ~0-98lto'o84 otherwise. When describing wde through a simple lin- 
ear function of the scale factor, our results do not evidence any time evolution. In the 
next few years new experiments will allow to measure the clustering of galaxies with 
a precision much higher than achievable today. Models like the one used here will be 
a valuable tool in order to achieve the full potentials of the observations and obtain 
unbiased constraints on the cosmological parameters. 

Key words: large-scale structure of Universe - cosmology: theory - cosmology: 
observations - cosmological parameters 



1 INTRODUCTION Universe is undergoing a phase of accelerated expansion 

driven by a new exotic component, dubbed d ark energy, 

The last decade was characterised by a dramatic change m , i •, ■ i , ^r^nv r .-, . . i HT- T^l 

whose ene rgy density is abou t TOyo or the total (IKiess et al.l 

our vision oi the Universe and by an impressive increase T^nr^a r, i ,, . , llir,r,rJ\ rm • i i ^i n ^ 

,,,,■' , , , r , 11998.: Pcrlmutter et al.lll999l l. This has been then conhrmed 

m size and quahty or available datasets. At the end or the i , , i . . i , , ■ • i i 

by other observati ons, such as the cosmic microwave back- 

twentieth century, the analysis or the luminosity-distance re- , ,„, ,„ [77^ — i 1 — nir,r,^nl [7: i — I — niorvonl 

. ■" ^ -^ ,,,,,, ground (CMB, e.g.,LH insha w et al.ll2003l : ISpergel et al.ll200a . 



lation in the Type la supernovae (SNIa) showed that the lS^r>rJ [TT , — ' — f~^ ^r,^^ r.,^T^ .11' 1 , . r 

■'^ ^ ^ ' 120071 : IKomatsu et al. 2009. 20111). the large scale structure of 

the U n iverse (LSS; e.g . lEfstathiou et al.ll2002l:|Percival et al.l 

l2002l : iTegmark et all |2004| : ISanchez et alT I2OO6I . bOOJ 
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iPercival et allboiol : iReid et alJlJOlol : iBlake et al.ll201ll ') and 
the numb er density of galaxy c lusters as function of their 
mass (e.g. IVikhlinin et all 12003 ). Different combinations of 
these probes have been used in the past years to constrain 
the dark energy equation of state parameter wde with about 
a 5-10% error. The analysis presented here aims at constrain- 
ing cosmological parameters and shed light on the nature of 
dark energy. 

Many present day and future galaxy redshift surveys, 
like the Baryonic Oscillation Spectroscopic Survey (BOSS , 
ISchlegel. White fc EisensteinI I2OO9I : lEisenstein et al.ll201ll ). 
the Panoramic Surv ey Telescope fc R apid Response Sys- 
tem (Pan-ST ARRS, [Ka iser ct all 120021 ') . the Dark Energy 
Survey (DES. lAbbottet al. 2005h. the Hob by Eberly Tele - 
scope Dark Energy Experiment fHE TDEX.lHiU et alll200i ') 
and the space based Euclid mission (|Laureiisll2009l V are de- 
signed to constrain the properties of dark energy, usually 
parametrised by its density, the present day value of w-q-e 
and its possible time evolution, with unprecedented preci- 
sion. This will help to reduce the number of possible models 
of dark energy that has been proposed in the last decade, 
that are either based on a time varying field (for a review 
see, e.g.. |Peebles fc Ratrall2003l ). or on modifications of the 
equation of general relati vity (modified gravity, for a review 
see, e.g.. iTsuiikawalbOlOl 'l. 

In this work we concentrate on the distribution of 
galax;ies on large scales as a means to constrain the his- 
tory and composition of the Universe, analysing it sta- 
tistically through the power spectrum, the Fourier trans- 
form of the two point correlation function. The shape 
of the galaxy power spectrum and correlation function 
contains information about the composition of the Uni- 
vers e and the non-linear evolution of clustering, bias 
(e.g. lMcDon"aidll2006l : iMatsubaralboOSbl : | jeong fc Komatsul 

Scoccimarrd 20041: 



( Crocce fc Scoccimarrd 20081: 
120081 : ISmith. Scoccimarro 



2009) and redshift space dist orti ons (e.g. 



spac ,^ 

Cabre fc Gaztafiagaj l2009al rbl: iJennings. Baugh fc Pascolil 
2OIII : iReid fc White! I2OIII '). Biasing is due to the fact that 



the objects that we observe do not trace perfectly the un- 
derlying dark matter distribution and redshift space dis- 
tortions are introduced when inferring the distance of a 
galaxy from the measured redshift, which is a sum of a 
cosmological component and a Doppler shift due to the 
peculiar motion of the emitter. The two-point statistics 
contain also a feature, called baryonic acoustic oscillations 
(BAO), that has been advocated as a powerful tool to probe 
the curvature of the Universe and that has been inten- 
sively studied over the past years. BAOs are the relic sig- 
nature in the matter distribution of the acoustic oscilla- 
tions in the baryon-photon plasma in the hot young Uni- 
verse. The BAOs show up in the correlation function as 
a quasi gaussian bump at scales r ~ 100 — 110/i^^Mpc 
iMatsubaral 12004 ) and as a sequence of damped quasi- 
harmonic oscillations at wave-numbers 0.01 feMpc ~^ < 

fc < 0.4/iMpc~^ in the power spectrum (|Sugivamal Il995l : 



lEisenstein fc Hu|[T998. 1999i). B AOs i n the CMB where pre- 
dicted by iPeebles fc Yul (|l970l ) and ISunvaev fc ZeldovichI 
(Il970h and first detec ted in the gala. x y dis tribution by 
ICole et all (|2005l ') and lEisenstein et all iJ2005l 'l. The BAO 
scale, as measured from the CMB, depends only on the 
plasma physics prior to recombination, which allows in 
principle to use it as a standard ruler. Non linear evolu- 
tion damps and shifts by few percent the acoustic peaks 



_.. iSanchez. Baugh fc Anguld 
Shethll2008l ). The full shape of 
the power spectrum and of the correlation function, as well 
as the BAOs alone, have been used, alone and in combina- 
tion with other indepe ndent datasets, to constrain cosmolog- 
ical parameters ( e.g.. IPercival et al.ll2007l: ISanchez fc Cole 



Cabre fc Gaztafiagaj l2009al: Gaztafiaga. Cabre fc Hui 
lcx„^u^„ ^+ „i llonnnl. It<'„„;„ ^<^ ^i Horn nl. l"D^,-^;,,.^i ^+ „1 



2003; 'Sanchcz ct al."2009'; 'xazin ot al."2010'; 'Porcival ct al 
[2010; Rcid ct aL,2010; Blake ct al. 2011; Tinker ct al. 201lJ) 
In order to obtain unbiased constraints from the infor- 
mation encoded in the large scale structure of the Universe, 
non linear distortions, bias and redshift space distortions 
need to be accounted for. Perturbation theory (PT, see 
iBernardeau et al.ll2002l . for a review) can successfully model 
the shape of the power spectrum and of the correlation 
function at redshifts larger than 2 = 1 or at large scales 
when includ ing at least third order term s in the density 
fluctuations (jjeong fc Komatsull2006l . l2009l ). In the past few 
years many groups have proposed different improvements 
over perturbation theory (e.g ., ICrocce fc Scoccimarro 



2006al l 
20071: 



JMcDonaldi 



Matarrese fc Pietron: 



Bernardeau. Crocce fc Scoccimarrd 
Matsubaral 



Matarrese fc Pietroni 



12008 



Pietronil l2008l: iTaruva fc Hiramatsul 

Smith. Hernandez-Monteagudo fc SeliakI 



12008 
l2008al ll 



Taruva et all l2009l : lEUa et all 120101). Wh ile 



this approaches, like for in s tance iMatsubaral (|2008d ) and 
iTaruva. Nishimichi fc Saitd (|2010h . can account for bias 
and redshift space distortions, most of them describe 
only the clustering of dark matter in real space. For this 
reason phenomenological approaches based on a given 
flavour of PT have been put forward in order to model the 
dark matter halo and galaxy distributions in real and in 
redsh i ft space (e.g., McDonald! 2006!: Crocce fc Scoccimarro 



20081 : ISanchez. Baugh fc Anguld I2OO8I: iJeong fc Komatst 
2009! : iMontesano. Sanchez fc Phleps! I2OIO! . the latter will 
be thereafter referred to as MIO). The model proposed 
by M IO is inspired by renormalis ed perturbation theory 
fRPT lCrocce fc Scoccimarrdl2006a! ld) and can describe the 
full shape of the mildly non-linear the dark matter and 
halo-power spectrum from numerical N-body simulation 
and allows to obtain unbiased constraints on the dark 
energy equation of state parameter. 

In this analysis we apply the model developed and 
tested in MIO to the power spectrum measured from the lu- 
minous red galaxies (LRGs), released publicly in the seventh 
data release of the Sloan Digital Sky Survey (SDSS DR7). 
Combining this measurement with the latest results from 
CMB, SNIa and the precise determination of the local Hub- 
ble constant. Ho, we obtain cosmological constraints for five 
different parameter spaces. The LRG sample and the mock 
catalogues used to estimate the covariance matrix of the 
data are presented in Section [2| The power spectra and co- 
variance matrix computed from the real and mock datasets 
are shown in Section [3| The CMB and SNIa datasets and 
the Ho measurement are described in Section 2) Section [5. II 
gives a short description of the model that we apply to the 
LRG power spectrum. In Sections 15.21 and 15.31 we illustrate 
the parameter spaces that we explore in our analysis and the 
technique used to perform the fits. In Section [5.4l we test our 
model against the mean power spectrum of the mock cata- 
logues and of the LRG, in order to determine whether it can 
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Table 1. Cosmological parameters and specifications of the 
LasDamas-Oriana simulations 
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Figure 1. Spline fit to the redshift distribution of the LRGs 
(dot-dashed line) and the original and modified n{z) of the mock 
catalogues (dashed and solid lines, respectively) 



describe accurately the full shape of the two-point statistic 
when a complex geometry is used. Our main results, the 
cosmological constraints for the five parameter spaces, are 
discussed in Section [6] and then compared with some of the 
most recent works on the topic in Section [T] Finally we sum- 
marise our results and draw our conclusions in Section [H] 



2 THE GALAXY SAMPLE AND THE MOCK 
CATALOGUES. 

In this section we describe the galaxy sample (Section 12. ip 
and the mock catalogues (Section 12.21) that we use in this 
work. 



2.1 



The luminous red galaxy sample from the 7th 
data release of SDSS 



Data release 7 (DR7, lAbazaiian et al.ll2009l ) is the last data 
release of the second phase of SDSS, know as SDSS-II. 
From the 929,555 galaxies, whose spectra have been mea- 
sured, we use the subs ample of lumi nous red galaxies (L RGs, 
lEisenstein et al.ll200ll ) presented in lKazin et al.l (|201CI ) and 
publicly availablqj. The catalogue contains 89,791 LRGs, in 
the redshift range 0.16 < z < 0.44 {z = 0.314), from the 
large contiguous area of 7646 deg'^ in the Northern Galactic 
Cap. The full survey also includes three equatorial stripes, 
that we do not consider. This causes a loss of less than 10% 
in galaxy number and volume, but the resulting geometry is 
simpler. Furthermore, the use of the Northern Galactic Cap 
only allows us to obtain a more accurate estimate of the 
statistical errors than for the full survey (see Section [2. 2p . 
The dot-dashed line in Figure [T] shows a smooth spline fit to 
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the sample number density as function of redshift. Together 
with the LRGs catalogue, we use a random one with about 
fifty times more objects, designed to reproduce the geom- 
etry and completeness of the galaxy sample and to have a 
radial number density proportional to the dot-dashed line 
in Figure [l] 



2.2 The moclt catalogues 

In order to test our analysis technique and to estimate the 
covariance matrix associated to the LRG power spectrum, 
we use the LasDamas mock catalogues (McBride et al., in 
prep.). The mocks have been constructed from a suite of 40 
large dark matter N-body simulations, dubbed Oriana, that 
reproduce a part of a universe characterised by a geometri- 
cally flat cosmology dominated by a cosmological constant 
and cold dark matter (ACDM). The cosmological parame- 
ters and specifications of the simulations are listed in Ta- 
ble [T] From each simulation a halo cata logue is extracted 
using a Friend-of-Friend algorithm (FoF. lDavis et aLlligSST l 
with linking length 0.2 times the mean inter-particle separa- 
tion. In order to match the LRG clustering signal, the haloes 
have then been populated with mock galaxies using a halo 



occupation distribution (HOD, iBerlind fc Weinberg 20021) 
within the halo model approach (HM, see lCoorav fc ShethI 
i2003 . for a review). The HOD parameters have been chosen 
in order to reproduce the galaxy number density and the 
projected correlation function of the observed SDSS DR7 
samples. From each simulation two (four) mock catalogues 
of the full SDSS DR7 volume (Northern Galactic Cap only) 
have been extracted. These mock catalogues, together with 
the mocks from two smaller companions of Oriana and the 
corresponding randoms, are publicly available Q 

In this work we use the 160 mock catalogues of the 
LRGs in the Northern Galactic Cap region. We modify the 
mocks and the corresponding random catalogue, which have 
the radial number density shown by the dashed line in Fig- 
ure [T] in order to reproduce the one of the LRG: the result- 
ing n{z) is indicated by the solid line. The mock catalogues 
contain on average 91137 galaxies. 
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3 THE POWER SPECTRA 

From the dataset and the mock catalogues just described we 
compute the power spectra and the covariance matrix that 
are used in the rest of the analysis. They are presented in 
this section. 



3.1 The LRG powrer spectrum 

To compute the power spectrum and the window function 
we need to convert the angular positions and redshifts of 
the galaxies and of the random points into comoving coor- 
dinates. This is done first inferring radial distances from the 
measured redshifts and then converting the spherical coor- 
dinates into cartesian ones. To do the first step we assume as 
fiducial the cosmology of the LasDamas simulations, shown 
in Table [B 

We compute the power spectrum, as well as the 
sur vey window function, using t he es timator introduced 
bv iFeldman. Kaiser fc Peacocki (|l994l . thereafter FKP). 
IPercival. Verde fc Peacockl ( 2004 . thereafter PVP) proposed 



a modification of the FKP approach to take into account the 
relative biases between populations with different luminosi- 
ties. In Appendix IBll we show that, thanks to the fact that 
the LRG sample is almost volume limited and composed by 
a relatively homogeneous class of galaxies, the shape of the 
power spectra recovered with the two methods are in ex- 
cellent agreement at linear and mildly non-linear scales. In 
Appendix|X]we summarise the most important equations of 
both estimators. 

At first we correct the gal axy catalogue for t h e loss of 
objec ts due to fibre collisions (jZehavi et al.l [2003; iMasiedil 
|200a). The SDSS spectrographs are fed by optical fibres 
plugged on plates, which forces the fibres to be separated 
by at least 55". It is then impossible, in a single expo- 
sure, to obtain spectra of galaxies nearer than this angu- 
lar distance. The problem is partially alleviated by mul- 
tiple exposures, but it is not possible to observe all the 
objects in crowded regions. Assuming that in a given re- 
gion of the n galaxies that satisfy the selection criteria we 
can measure only m ^ n redshifts due to fibre collision 
and assuming that the missed galaxies have the same red- 
shift distribution of the observed ones, we assign to the 
latter a weight Wi — n/m. This ensures that the sum of 
the weights in a given region of the sky is equal to the 
number of selected galaxies n. Secondly to each LRG and 
random object at position x, where the number density is 
n(x), we associate a weight w{x) = (1 -|-p„n(x))^^, with 
Pw = 40000/i~''Mpc'^. This value has been chosen in order 
to minimise the variance of the measured power spectrum 
in the range 0.02 /iMpc"^ < A; < 0.2/iMpc"\ In Appendix 
IBll we show the results of the tests to analyse the impact of 
different choices of pw and corrections, namely fibre collision 
and completeness, on the recovered power spectrum. 

To compute the power spectrum we assign the LRGs 
and the random objects, weighted as described before, to a 
cubic grid with N = 1024^ cells and side L = 2200 /i"^ Mpc 
using triangular shaped cloud (TSC) as mass assignment 
scheme (MAS). For each cell, we compute the F(pi) field of 
equation HAla|l . We then perform the fast Fourier transform 
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Figure 2. LRGs power spectrum (blue dots connected with solid 
line) and corresponding l-a error bars from the mock catalogues 
(shaded area). The green dashed and the red solid lines show, 
respectively, the linear and model power spectra computed using 
the mean value of the cosmological parameters of the ACDM 
cosmology show in the last column of Table |3] 



(FFT) using the publicly availa ble software fft'VvIj (Fa stest 
Fourier Transform in the West. lFrigo fc Johnsonll2005l ). We 
correct each Fourier mode by dividing it by "^^ \W(k — 
2fcNn)p, where W^(k) is the Fourier transform of TSC, /cn — 
■k{N^^^ /L) is the Nyquist wavenumber and n is a 3D integer 
vector. Finally we spherically average the Fourier modes and 
subtract the shot noise (equation I A5a|) . 

The window function is evaluated similarly. We assign 
the objects of the random catalogue to four cubic grids 
with N = 1024^ cells and sides L = 2200, 4400, 8800 
and 17600 h~^ Mpc and compute the field G'(x) of equation 
(|A6a[l . We use the four grids with different dimensions in 
order to be able to compute the window function up to very 
large scales. For each box, we perform the FFT, correct the 
Fourier modes, spherical average and subtract the shot noise 
(eguation lATap . For each window function G^{k), we discard 
all the modes with wave-number k > 0.65 ^n and, when two 
or more window functions overlap, we consider only the one 
computed in the larger volume. This choice is motivated by 
the fact that, for a given band in wavenumber, the larger 
volume window function has been computed averaging over 
a larger number of modes than the ones from smaller vol- 
umes. Finally we merge the four window functions in order 
to obtain a single curve. 

According to equation HA4|) . the observed power spec- 
trum Po(fc) just described is a convolution of the "true" 
power spectrum Pt{k) with the window function. This con- 
volution is computationally time consuming, in particular 
when it has to be performed repetitively. Therefore we trans- 
form this convolution into a matrix multiplication: 



Po(fcO = ^ W^(fci,fc,)Pt(fc,) - CG'(fcO 



(1) 
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Figure 3. Panel a): power spec tra of Figure [2] divided by a linear 
power spectrum without BAOs IjEisenstein &: Hulll998r ). Panel b): 
rows of the window matrix corresponding to the fc-bands of the 
measured LRG power spectrum. For clarity only one every third 
row is shown. 
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Figure 4. Panel a): mean power spectrum (blue dots connected 
by solid line) with \-a variance (blue shaded area) from the mock 
catalogues. The linear and model power spectra, convolved with 
the window function, are shown with green dashed and red solid 
lines, respectively. Panel h): same po wer spectra divided b y a 
linear power spectrum without BAOs JEisenstein fc Hulll998l) . 



W{]ti,kj) = ajfc|/j^^dcos(6')G^(|ki — kj|) is the window 
matrix normalised such that J]] Wiki, kj) = \ \li. The coef- 
ficients aj corresponding to the wavenumber fcj are derived 
using the Gauss-Legendre decomposition. The second term 
in the right hand si de arises from the integral constraint 
(jPercival et al.ll2007n where C is a constant determined by 
requiring that Po(0) = 0. 

The LRG power spectrum is shown in Figure [2] with 
blue dots connected by a solid line. We also show the lin- 
ear (green dashed line) and the model (red solid line) power 
spectra computed from the best fit parameters obtained as- 
suming a flat ACDM cosmology (see section [6711 for more de- 
tails) . Both of the power spectra have been convolved with 
the window function as in equation ([ij and their amplitude 
has been boosted in order to match the observed one. The 
blue shaded area represent the variance from the mock cata- 
logues. Panel a) of Figure |3] shows the same quantities, with 



the same colour and line coding, but divided by a smooth lin - 
ear power spectrum without BAOs (jEisenstein fc HulllQQa '). 
Panel b) of Figure [3] shows one every third row of the 
window matrix Wiki^kj). Because of the relatively simple 
and uniform geometry of the sample used, the window func- 
tion has its maximum at fc = 0/iMpc~^ and decreases very 
steeply. This translates into the very sharp peaks at kj ~ ki 
in the window matrix rows, as shown in the figure. 



3.2 The mock power spectra and covariance 
matrix 

We compute the power spectra from the 160 realisations 
and the window function as we do for the LRG sample in 
the previous section. We then compute the mean P{k), the 
standard deviation and the covariance matrix C, whose ele- 
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Figure 5. Correlation matrix computed from tiie LasDamas 
mock catalogues. 



ments are defined by: 

Crur. = J-, — T Y. (Plikm) " P (kn^)) {Pi (kn) - P(fc„)), 

(2) 
where Piikm) corresponds to the measurement of the power 
spectrum at the m-th fc-bin in the l-ih reahsation. 

Panel a) of Figure [4] shows the mean power spectrum 
and its variance (blue dots with solid line and shaded area) 
obtained from the mock catalogues. The green dashed and 
the red solid lines show the linear and the best fit model 
power spectra (see Section [5.4. II for more details) convolved 
with the window function of the mocks and rescaled by the 
bias. Panel b) of Figure [4] shows the same power spectra of 
panel a) divided by a power spectrum without oscillations 

The correlation matrix of the mock catalogues, defined 
as Cmn/\/CmmCnn, IS showu lu Figure [5] The mode cor- 
relation caused by the convolution with the window func- 
tion is visible in particular near the diagonal. Non linear 
mode coupling is present at small scales and its strength in- 
creases with increasing k. Although the correlation becomes 
important for k > 0.2 h Mp c~^, it is not negligible already at 
k ~ .1 /iMpc~^. Recentlv lSamushia. Percival fc Raccanellil 
(|201ll ) showed that different methods of constructing the 
random catalogues can affect the estimated covariance ma- 
trix. 

The version of the mocks available when writing this 
article does not contain information about the luminosity 
of the galaxies, completeness and fibre collision. We cannot 
therefore test the impact of different estimators and cor- 
rections on the errors. However we analyse the impact of 
different pw on the power spectra and errors as measured 
from the mocks: the results are reported in Appendix IB2I 



4 ADDITIONAL EXPERIMENTS 

Later in this article we use the measurement of the LRG 
power spectrum described in the previous section to con- 
strain cosmological parameters. In order to obtain tight con- 
straints it is not possible to use this measurements alone, 
since some parameters have a weak dependence on the shape 



of the power spectrum and others present strong degenera- 
cies. Using the information coming from independent experi- 
ments, like the cosmic microwave background (Section 14. ip 
and the type la supernovae (Section l4.2|l , or prior knowledge 
of some of the parameters, like the local Hubble parameter 
(Section 14. 3p . it is possible to greatly improve the accuracy 
of the analysis. 



4.1 Cosmic microwave background 

Accurate measurements of the CMB temperature and po- 
larisation anisotropies provide a powerful tool to constrain 
cosmological parameters. In this work we use the observa- 
tions from five different instruments. 

The Wilkinson Microwave Anisotropy Probe (WMAP) 
satellite produced full sky maps of the CMB with a res- 
olution of 0.2°. We make use of the temperature angu- 
lar power spectrum in the multipole range 2^1^ 1000 
and the temperature-E polarisation cross power spectrum 
in the range 2 ^ I gC 450 from t he 7th year data re- 
lease (WMAP7. IJarosik et all I2OI0I : iKomatsu et~aLl I2OIII : 
iLarson et aLlbOllf T 

We also use measurements of higher multipoles from 
other four experiments which observe the CMB tempera- 
ture anisotropy on small patches of the sky with a much 
higher resolution. In order to avoid complex correlations, 
we use these experiments only for multipoles that do not 
overlap with the WMAP measurements. We use the mea- 
surements of the temperature angular power spectrum i) in 
14 bandpowers in the range 910 ^ / ^ 1850 fro m the Ar- 
cminute Cosmology Bolom eter Array (ACBAR, IKuo et al.l 
l2007l : lReichardt et al.ll2009l ). ii) in 6 bandpowers in the range 
855 ^ I ^ 1700 from the Cosmic Backg round Imager (CBI, 
for latest results see lSievers et al.ll2009l ). iii) in 7 bandpow- 
ers in the range 925 ^ Z ^ 1400 from the 2003 flight of 
the Balloon Observations Of Millimetric Extragalact i c Ra- 
diation and Geophysics (BOOME RanC '.lonos ct a ll I2OO6I: 
MacTavish ct al. 2006; Montrov et al. 2006; Piaccnti ni et all 
20061 ) and iv) in 11 bandpowers in the range 974 < I < 1864 



from the flnal da ta release of QUES T at DASI (QUAD, for 
latest results see lBrown et al.ll2009h . We also use E and B 
polarisation (EE and BB) and the cross temperature-E po- 
larisation (TE) angular power spectra measurements from 
the latter three experiment, i) CBI: EE in 7 bandpowers 
in the range 860 ^ I ^ 1800, BB in 5 bandpowers in the 
range ^ / ^ 5000 and TE in 8 bandpowers in the range 
860 < Z ^ 1800; u) BOOMERanG: EE and BB in 3 band- 
powers in the range 600 < I < 1000 and TE in 6 bandpowers 
in the range 450 < Z < 950; in) QUAD: EE, BB and TE in 
17 bandpowers in the range 488 < I < 1864. 



4.2 Type la supernovae 

Type la supernovae (SNIa) provided the flrst evidences of 
an accelerating Universe and the need of a new exotic com- 
ponent, called dark en ergy, to explain it (jRiess et al.lll998l : 
IPerlmutter et al.lll999l ). From the light curve, i.e. the lumi- 
nosity variation as function of time, of a SNIa it is pos- 
sible to measure the absolute luminosity, from which the 
distance to each object is inferred, probing in this way 
the distance-redshift relation. The first step is performed 
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thanks to models which encode intrinsic variations due to 
the physics of the SNIa, the effects from galactic and inter- 
galactic medium and selection effects in different ways. Be- 
cause of this, models produce different results, which impact 
the accuracy at which cosmological parameter s can be mea- 
sured and can introduce s ystematic effects (e.g. lHicken et al.1 
I2OO9I : iKessler et aLll2009l ) . 

SNfa data have been collected by many surveys de- 
signed according to different strategies and carried out us- 
ing a large variety of telescopes. Each of them observes a 
relatively small number of events, typically less than a hun- 
dred. In order to increase the number of objects, recently 
collections of SNIa from different surveys have been cre- 
ated. In this work we u se one of such samples, the Union2 
IjAmanullah et al.ll201Cl V It consists of 557 SN drawn from 
17 datasets in the redshift range 0.015 ^ z ^ 1.4 and is the 
largest available sup ernovae sample to dat e. It extends and 
improv es the Union fKowa lski et al.l 120081 ) and the Consti- 
tution (JHicken et al . 200!^) datasets: all the ligh t curves of 
the s elected SNIa have been fitted with salt2 (|Guv et al.l 



rameter spaces explored in Section[6]and the method used to 
extract cosmological information. Finally, we test the model 
against the LasDamas and the LRG power spectra (Section 
[53J. 



20071) and an improved analysis of systematics is presented. 
J (|2010f l showed that the inclusion of sys- 



AmanuUah et al.l 



tematic errors when fitting cosmological parameters, using 
only SNIa or in combination with independent probes (BAO, 
CMB and Hq), increases the associated errors leaving the 
best fit value almost unchanged. In Sections 16. II and 16. 31 we 
will further comment on the impact of systematics and of 
light curve fitters. 



4.3 Hubble parameter 

The Sup ernovae and HO fo r the Equation of State Program 
f SHOES. iRiess et al.ll2009^ aims at the direct measurement 
of the Hubble parameter at present epoch Hq to better 
than 5% accuracy. The SHOES team identified 6 nearby 
spectroscopically typical SNIa, that have been observed be- 
fore maximum luminosity, that reside in galaxies containing 
Cepheids and that are subjects to low reddening. Thanks 
to 260 Cepheids observed with the Near-Infrared Camera 
and Muhi-Object Spectrometer (NICMOS) on the Hubble 
Space Telescope (HST) in the 6 host galaxies and in the 
"maser galaxy" NGC 4258, the authors could calibrate di- 
rectly the peak luminosity of the SNIa. Combining these 6 
objects with 240 SNIa at redshift 2 < 0.1, they measure 
the Hubble parameter to be Ho = 74.2 ± 3.6kms~^Mpc^^. 
The error includes both statistical and systematic uncertain- 
ties. In this work we use this result as a prior knowledge of 
Ho, under the assumption that the associated likelihood is 
a gaussian with the mean and standard deviation measured 

by the SHO ES team. 

Recentlv lMoresco et al.l (|2010l ) measured, using spectral 
properties of early-type galaxies, the Hubble parameter to 
be Ho = 72.6 ± 2.9kms~^ Mpc~^ at 68% co nfidence level. 
While finishing this work, iRiess et all (|201lh refined their 
analysis obtaining Ho = 73.8 ± 2.4kms~^Mpc~^. 



5 METHODOLOGY 

In this section we describe the model for the full shape of the 
power spectrum that we use throughout the rest of this arti- 
cle (Section 15. ip . In Sections 15.21 and 15.31 we present the pa- 



5.1 The model 



MIO introduced a model for the full shape of the power spec- 
trum inspired by renormalize d perturbation theory (RPT, 
ICrocce fc Scoccimarrd l2006al | b|) and analogous to the one 



for t he correlation function of Crocce fc Scoccimarrd (|2008f ) 



and ISanchez. Baugh fc Anguld ( 20081 1. In this model they 
parametrize the non-linear power spectrum as 

P(fc,2) =b' (e-('='"'*)'PL(fc,2)H-4McPiioop(fc,2)) , (3) 

where the linear bias b, the damping scale of the BAO oscil- 
lations fc* and amplitude of the mode coupling contribution 
Amc are free parameters. Pi^{k,z) is the linear power spec- 
trum and 

Piloop(fc) = -^ Jd'q\F2{-k - q, q)i'PL(lk - ql)PL(g) (4) 

is the lowest order term arising from the coupling of two 
initial modes. The exponential damping and Piioop(fc) are 
the approximations of the non-linear propagator and of the 
mode coupling power spectrum described by RPT. 

The model of equation (|3J has been successfully tested 
in MIO against a suite of very large, median resolu- 
tion N-body nume r ical simulations, called L-BASI CC II 
l|Angulo et al.ll2008l : [sinchez. Baugh fc Anguld I2OO8I ). MIO 
computed power spectra for the dark matter distribution 
and several halo samples, selected according to their mass, 
both in real and redshift space at z=0, 0.5 and 1. They have 
shown that in all cases equation ((3)1 describes accurately 
these power spectra for fc < 0.15/iMpc~^ and that it allows 
to obtain unbiased constraints on cosmological parameters. 



5.2 Parameter spaces 

In the description of the datasets presented in Sections [2]|4j 
a strong assumption was implicitly made: the Universe is, 
at large scales, statistically homogeneous and isotropic and 
that density and velocity fluctuations around their mean 
values are small. In the following, we further assume that 
the fluctuations set by the initial conditions were adia- 
batic, gaussian and almost scale invariant and that they do 
not present tensor modes. The WMAP7 data, both alone 
and in combination with BAO, Ho and SNIa, conflrm that 
those assumptions ar e correct at the 95 % confidence level 
(Komatsu ct al. 2oTll: lEarson et aLlbOllI ). 

Within this framework we analyse five different cos- 
mological models, explored using five combinations of 
experiments: CMB, CMB-l-PLRG(fc), CMB-|-PLRG(fc)+-H'o, 
CMB-hPLRG(fc)+SNIa and CMB-hPLRG(fc)+-ffo+SNIa. 

The "concordance" ACDM model is the simplest model 
able to successfully describe a large variety of cosmological 
datasets. It describes a geometrically flat (fik ~ 0) uni- 
verse with a cosmological constant A, whose equation of 
state parameter uide = — 1 is constant in space and time, 
and pressureless cold dark matter (CDM) as main compo- 
nents. This cosmology can be characterised by six parame- 
ters: the baryon and dark matter densityies ojb = f2b/i^ and 
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cjdm ~ 0.mAh? , the scalar spectral index Us and the ampli- 
tude Ab of the primordial fluctuations, the optical dept r to 
the last scattering surface, assuming instantaneous reionisa- 
tion, and the ratio between the horizon scale at decoupling 
and the angular diameter distance to the corresponding red- 
shift 9: 

^ACDM = {ijOi,,ujum,ns,Aa,T,Q\h,ki,,Amc,Asz)- (5) 

The four parameters after the semicolon are related to the 
modelling of the matter power spectrum (&, ki, and Amc 
from equation [3]) and of the CMB angular power spec- 
trum {Asz, amplitude of the contribution to the CMB at 
large I from the Sunyaev-Zeldovich effect). These parame- 
ters are marginalised over when showing the cosmological 
constraints in the section [S] 

If we then drop one or both assumptions on geometry 
and on the value of the dark energy equation of state, we 
obtain three cosmologies characterised by the following pa- 
rameter spaces: 

• variable curvature, uide = —1: 

^fcACDM = (a;b, a;DM, fik, Ws, A, r, 9; 6, A;*, ^mc, ^Sz); (6) 

• zero curvature, tode = const: 

SwCDM = (cJb,^DM, WDE,ns,y4s,r, 6;b, A:*, Amc^sz); (7) 

• variable curvature, wde = const: 

^kwCDM = (l^bi^^JDM, r^kjlt'DEi^.s, ^s, T, 6; b, fc,, AmC, Asz)- 

(8) 

As last case, we consider a flat Universe in which «;DE(a) 
ev olves with time. We adopt th e pa rametri s ation proposed 
bv lChevallier fc Polarskil (|200ll ) and lLindeJ (|2003l '): 



WDE(a) — Wo + TOa(l — O) . 



(9) 



Although not physically motivated, it can describe accu- 
rately a big variety of equations of state derived from scalar 
fields with the use of only two parameters: its value today. 
Wo, and its first derivative with respect to a, tUa- The result- 
ing parameter space is: 

SwaCDM = {u}h,^oM,wo,W!i,ns,AB,r,Q;b, k-^,, Amc, Asz)- 

(10) 
Other cosmological quantities can be derived from the 
ones just listed. In particular we are interested in: 



^dcr = (Ode, f^M, -ff, i0,O"8, Zrc). 



(11) 



The density of dark energy, f^DE, is obtained from a com- 
bination of r2k, i^M = oJb + i^DM and B. From there, 
the total matter density Om = 1 — Slk — f^DE, the Hub- 
ble parameter h = ^ujm/^m and the age of the uni- 
verse to = /„ da/{aH{a)) are derived. The present day 
rms of linear density fluctuation in a sphere of radius 
8/i~^Mpc, as, is computed from Ab- From r, H, fib and 
SlpM it is possible to estimate the redshift of reionisation Zrc 
(JTegmark. Silk fc Blanchardll 19941 ). 

5.3 Practical issues 

To constrain the sets of cosmological and model parameters 
just described, we use the Markov Chain Monte Carlo 
technique (MCMC, ICilks. Richardson fc Spieeelhalteil 



Table 2. Prior ranges for the primary cosmological and the model 
parameters 



Parameter 



lower limit upper limit 



un, 



: Uuh'^ 



f^k 

'i'DE ('"'o) 

Wa 

Jla 

log(lOio^s) 

r 

1000 



0.005 


0.1 


0.01 


0.99 


-0.3 


0.3 


-2 





-2 


2 


0..5 


1.5 


2.7 


4 


0.01 


0.8 


0.5 


10 



fc* 


0.01 


0.35 


^MC 





5 


Asz 





2 



119961 : IChristensen fc Meveil l2000l ) as implemented in 
the free software C OSMOMqj (Cosmological MonteCarlo, 
iLewis fc Bridle! I2OO2I '). The CMB and linear matter power 
spectra are computed with a modified version of CAMB 
(Code for Anisotropics in the Microwave Background, 
iLewis. Challinor fc L ascnbv 200(r ) that allows to consider 
time varying wde |j(Fang. tlu fc Lewis! !2008l ). For each 
choice of parameter space and probes we run eight in- 
dependent chains . The ir execution is stopped when the 
!Gelman fc Rubin! ^199^ 1 criterion R < 1.02 is satisfied. The 
MCMC requires some prior knowledge of the parameter 
space that is explored. We assume for all the primary 
parameters (equations I5I8I and I10|l flat priors in the ranges 
listed in Table (2] The model parameter b is analytically 
marginalised over a n infinite flat prior (equation F2 in 
!Lewis fc Bridle!!2002! ). 

All likelihoods used to compare the results from the 
cosmological probes described in sections l2.1ll4.1ll4.2l and l4.3l 
with the corresponding models are assumed to be Gaussian. 
In the case of the LRGs power spectrum, this is 



£cxexp(--x^(e) 



where 



x'W 



(d- 



■ mfc- 



\d-m) 



(12) 



(13) 



is the standard x^, in which d is an array containing the 
band power P{k), t{9) contains the model computed for the 
set of parameters 9 and then convolved with the window 
function and C is the covariance matrix of the measurement 
presented in section [32] 

In section 12.11 we assume a fiducial cosmology in order 
to convert redshifts and angles to physical coordinates. Dif- 
ferent choices of these parameters result in modifications of 
the measured LRG power spectrum. The ideal case would be 
to recompute, for each step of the MCMC chain, the power 
spectrum and the window function according to the given 
cosmology, but this is not computationally feasible. Under 
the assumption that the survey covers a wide solid angle. 



http : //cosmologist . inf o/cosmomc/ 1 



http : //camb ■ Inf o/ppf TJ" 
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it is possible to incorporate these distortions in the correla- 
tion function when computed using differe nt cosmologies by 
rescaling the distance scale by the factor (JEisenstein et alj 
l200d ): 



Dl 



\zu 



D^^i^n 



(14) 



where the effective distance Dv{zm) to the mean redshift of 
the sample Zm is, 



Dv{z) = 



-Di(Zni) 



H{Zm) 



1/3 



(15) 



with DA(zm) the comoving angular diameter distance. Since 
the power spectrum is the Fourier transform of the correla- 
tion function, this holds also for the former. Thus, at each 
step of the chain we multiply the wave-number of the model 
power by a and its amplitude by 1/a^ in order to rescale it 
to the fiducial cosmology. 

5.4 Testing the model 

In the first part of this section we extend the analysis of 
MIO and fit our model against the mock catalogues in or- 
der to test wether it provides an accurate description of the 
measured power spectrum and unbiased constraints of the 
dark energy equation of state parameter also when a com- 
plex geometry is involved. We also test the stability of the 
cosmological parameters in the wCDM case as the maxi- 
mum k included in the analysis varies, when CMB and the 
PhB.a{k) are combined. 



5.4-1 Mock catalogues 

In this section we follow MIO and assume all the parameters 
fixed, except for tode. Under this assumption, equation (|14|) 
links univocally variations of the dark energy equation of 
states to stretches a of the model power spectr um. There- 
fore w e consider the latter as a a free parameter (jHuff et al.l 
[2003). 

Using a MCMC approach, we explore the parameter 
space defined by ^ = (fc*, 6, Amc, q)- We chose priors with 
a constant probability within the following ranges: 

• 0/iMpc~^< fc* < 0.35/1 Mpc"\ 

• s; Amc < 10, 

• 0.5 s; Q < 1.5. 

The bias b is marginalised analytically over an infinite flat 
prior as described in section [531 

Figure [6] shows the one-dimensional marginalised con- 
straints on the parameters fc*, Amc and a varying the maxi- 
mum value of the wave number fcmax of the measured power 
spectrum used to perform the fit; we keep the minimum k 
fixed to 0.02 /iMpc~^. The blue circles connected by solid 
lines and the shaded areas correspond to the mean and the 
standard deviation of these parameters. For every fcmax the 
model is evaluated for k ^ 0.2/iMpc~^ and then convolved 
with the window function. Since each row of W{ki,kj) is 
sharply peaked at kj ~ ki, the main contribution to Po(fci) 
comes from modes near ki: therefore the constraints shown 
in Figure |6] depend weakly on the exact wave number range 
in which the model is computed, as long as it is larger than 
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Figure 6. One-dimensional marginalised constraints on the pa- 
rameters fc*, j4mc and a as function of the maximum value of 
ki (feniax) used to fit the model of equation ^ to the LasDamas 
mean power spectrum. The mean values and the standard devi- 
ation are indicated, respectively, by circles connected with solid 
lines and shaded areas. The maximum value of k at which the 
model is computed before the convolution with the window func- 
tion is fc = 0.2/iMpc-i. 



the range of the measured power spectrum that we fit. As 
fcmax increases, more modes are included in the fit and the 
errors decrease. The constraints on a are compatible with 
unity for all the fcmax considered. Considering that the vol- 
ume of the simulations used in MIO is much larger that the 
one sampled by the LRGs and consequently that the errors 
in the former case are smaller than in the latter, we decide 
to further consider scales 0.02/iMpc~^ ^ fe ^ 0.15/iMpc~^. 
The constraints on fc* and Amc exhibit, respectively, a 
monotonic increase and decrease. As explained in MIO, the 
approximate mode coupling power in equation ([Sjl is about 
30% larger than the exact value at fc ~ 0.15 — 0.2/iMpc~^: 
this forces Amc to decrease to ~ 0.7 — 0.6 and A:* to in- 
crease in order to maintain the shape of the resulting power 
spectrum unvaried. 

Figure [7] shows the two-dimensional marginalised con- 
straints in the fc* — a plane as obtained from the mock cat- 
alogues for fcmax = 0.15/iMpc~^. The inner dark and outer 
light shaded areas represent regions whose volumes are 68% 
and 95% of the total likelihood. This representation of the 
two-dimensional constraints will be used through all section 
[6] The independence of a from fc* or Amc makes the con- 
straints on the former robust. The latter two parameters, 
instead, are strongly degenerate, as it is possible to describe 
accurately the overall shape of the power spectrum compen- 
sating an increase (decrease) of fc* with a decrease (increase) 
of Amc- 

The model power spectrum indicated by a solid line in 
the two panes of Figures 3] has been computed using the 
best fit values of the parameters as obtained in this section 
for fc,„ax = 0.15 /iMpc-^: fc* = 0.26/iMpc-\ Amc = 0.61. 
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1.05 



e 1.00 



0.95 




0.20 
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0.25 
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Figure 7. Two-dimensional marginalised constraints in the fe* — a 
plane obtained applying the model of equation ((Sjl to the Las- 
Damas mean power spectrum. The results are show for fcmax = 
0.15/iMpc^^. The inner dark and outer light areas represent the 
68% and 95% confidence level, respectively. 



The bias has been computed maximising the likelihood of 
equation (|12p with all the other parameters fixed. 



5. 4 -2 The Luminous red galaxy sample 

After testing the robustness of our model at mildly non- 
linear scales, we test here the dependence of the cosmolog- 
ical parameters upon fcmax. We use the wCDM cosmology 
defined in equation (O and the combination of the LRG 
power spectrum and CMB measurements. Figure |S] shows 
the one-dimensional marginalised constraints on the param- 
eters oJb, ttJDM, wde, f^DE, ^M, o"8 and h as function of fcmax- 
The circles connected by solid lines and the shaded areas 
show the mean and the standard deviation as obtained from 
the MCMC. The model is computed for k < 0.2/iMpc"^ 
and then convolved with the window function; as before, 
the measured parameters are mostly insensitive to this limit. 
Although some parameters are more stable with respect to 
changes of fcmax than others, there are no significant trends 
or deviations. The error-bars of the parameters decrease as 
fcmax increases: we interpret this as a sign of new information 
coming from these scales. 

We also test the impact of using CMB measurements 
from WMAP only. We obtain the same constraints, but with 
errors larger by 5-10%, due to the loss of information at 
small angular scales, i.e. large multipoles I. Appendix IB3I 
describes how the results just described depend upon pw 
and Wi. A more extensive analysis of the wCDM cosmology 
is presented in section [ 
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Figure 8. One-dimensional marginalised constraints of the 
wCDM parameter space on the parameters wj,, tJuMi '^DEi ^DEi 
Cm, erg and h as function of the maximum value of k (fcmax) 
as obtained when combining the LRGs power spectrum with the 
CMB data. Circles connected by solid line and shaded areas are, 
respectively, the mean and the 68% confidence level. The model 
is computed for fc < 0.2hMpc~^ and then convolved with the 
window function. 



6 RESULTS: THE COSMOLOGICAL 
PARAMETERS 

In this section we present the constraints on the cosmologi- 
cal parameters obtained from the five cosmological scenarios 
introduced in section [5^ 



6.1 The concordance cosmology 

The flat ACDM cosmology, parametrized by the first six 
quantities of equation ([5]), is the minimal model able to de- 
scribe data coming from many independent probes. Table [3] 
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Table 3. Marginalised constraints on the cosmological parameters of the ACDM parameter space from the combination of 
probes listed in the header of the table. The quoted values are the means and widths of the posterior distribution containing 
68% of the total area. 
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Figure 9. Two-dimensional marginalised constraints of the 
ACDM parameter space in the Odm — Ho plane. Blue, red and 
green shaded areas enclosed in dashed, solid and dot-dashed 
lines shows the constraints from CMB alone, CMB+PlrgC^) and 
CMB+PLRG(fc)+Ho+SNIa, respectively. The inner darker and 
the outer lighter areas are the 68% and 95% confidence level. 
The black dotted line that runs across the figure shows the locus 
defined by equation QmHq ~ const. 



summarises the complete list of one-dimensional constraints 
for the primary and derived cosmological parameters ob- 
tained in this section. We quote the mean values and width 
of the posterior distribution containing 68% of the total area, 
which for a Gaussian distribution corresponds to the stan- 
dard deviation, as obtained for the five different combina- 
tions of the four experiments used in this work. This con- 
vention will be followed in the rest of this article. 



The CMB experiments described in section \4A\ provide 
measurements of the temperature and polarization angular 
power spectra with very high accuracy. The blue shaded ar- 
eas enclosed in dashed lines in Figure [9] show the 68% and 
95% confidence level in the JIdm— J^o as obtained when CMB 
information only is used. The apparent position of the peaks 
in the CMB power spectra is proportional to their physical 
scale, which depends on the composition of the early uni- 
verse (baryons, dark matter and radiation), and the angular 
diameter distance to the last scattering surface, a function 
of Ho and of the density and equation of state parameter of 
matter, dark energy and of curvature. Since here we consider 
a flat geometry and fix woe to -1, a degeneracy between 
the matter densit y and the Hubble pa rameter appears. It 
has been shown bv lPercival et al.l (120021 1 that this effect, to- 
gether with the preservation of the relative amplitude of the 
peaks, leads, in a ACDM universe, to a degeneracy along 
the curve defined by r^M^i^ — const, which is displayed in 
Figure |9] by the dotted line. The accurate detection of the 
third peak in the temperature power spectrum, whose rela- 
tive amplitude with respect to the first two is proportional 
to the matter-radiation ratio, helps reducing this degener- 
acy. In this case we measure Hm = 0.259 ± 0.026, Ho = 
71.7tl:t kms"i Mpc-i and ujom = (10.96 ± 0.52) x 10"^ 

The inclusion of information from the large scale struc- 
ture can break or reduce some of the degeneracies in the 
CMB. The shape of the power sp ectrum depends upon 
i^Mh and, more weakly , Slb/fiM (e.g. jEfstathiou et al.ll2002l : 
ISanchez fc Cole Il2008h . Thanks to this, the errors on the 
cosmological parameters decrease up to about 30%. In par- 
ticular, we measure a decrease in the allowed region by 
about one third in the three parameters considered be- 
fore: ^M = 0.27 ± 0.018, Ho = 70.8 ± 1.6 kms^^Mpc-i 
and oj-DM ~ (11.23^0 35) ^ 10^^. The two-dimensional con- 
straints in the Om — Ho plane are shown in Figure [5] with 
the red shaded areas within solid lines. 



When the SNIa and Ho measurements are also used, 
we obtain ^m = 0.265tJ5;j;i^, Ho = 71.2 ± 1.4kms-i Mpc'^ 
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Figure 10. Two-dimensional marginalised constraints of the 
ACDM parameter space in the f2m — erg plane. The colour and 
line coding is the same as in Figure |9] 



Table 4. Marginalised constraints on the cosmological pa- 
rameters of the ACDM parameter space from the combi- 
nation of CMB-|-PLRG(fc)+^0+SNIa when systematic er- 
rors are not considered and when the mlcs2k2 SNIa light 
curve fitter is used. 
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and CJDM = (ll-13la32) x 10~^, which means a 10-15% 
increase in accuracy. The two dimensional 68% and 95% 
confidence levels for the former two parameters, when all 
the four probes are used, are shown in Figure [9] by the green 
shaded areas enclosed by dot-dashed lines. 

Figure [To] shows the two-dimensional marginalised con- 
straints in the SIm — os plane for the same combination of 
datasets; colour and line coding is the same as in Figure 
[9] The correlation between the two parameters is caused 
by the fact that an increase (decrease) of Q.m causes a de- 
crease (increase) in the amplitude of the power spectrum 
that can be compensated by a larger (smaller) value of ag. 
For the three cases shown in the figure we find that the 
one-dimensional constraints on the latter parameter are, re- 
spectively, (78 = 0.796 ± 0.027, 0-8 = 0.807 ± 0.022 and 
as = 0.802 ± 0.021; the error on the latter two decreases 
by 20 and 22% respectively to the CMB only result. 

The solid lines in Figures [2] and [3] show the model 
power spectrum computed using the cosmological param- 
eters listed in the last column of Table [3] and the best fit 
model parameters fc* = 0.28/iMpc~^ and Ayic ~ 0.72, as 
obtained from the MCMC for fcmax = 0.15/iMpc"\ The 
bias has been computed maximising the likelihood of equa- 
tion (|12|l with all the other parameters fixed. 



pie consists of 288 supernovae from 5 different experiments 
whose light curves have been fitted both with the mlcs2k2 
(|jha. Riess fc Kirshneill2007l l and SALt2 models. Since most 
of the objects in SDSS SN are also part of the Union2, we 
do not expect large differences between the two samples if 
SALT2 is used. Therefore we show here only the results ob- 
tained using the mlcs2k2 fitter. The one-dimensional con- 
straints that we extract from the combination of all four 
probes, when the systematics in the Union2 are ignored and 
when the SDSS SN with mlcs2k2 are used, are listed in Ta- 



ble]^ As expected, and in agreement with I Amanullah et al.l 
(|2010f ) , ignoring systematic errors globally reduces the errors 
on the recovered values of the parameters without chang- 
ing sensibly the mean value. As example we obtain Om ~ 
0.266 ±0.013 (14% decrease with respect to the correspond- 
ing case in Table [SI, Ho = 71 ± 1.2kms~^ Mpc~^ (12%), 
OJDM = (I1.16lj;;29) X 10"^ (14%) and ag = 0.803 ± 0.02 
(5%). On the other hand, the use of mlcs2k2 changes the 
posterior distribution sensibly, without infiuencing its width. 
We measure that Om, <j-'dm and Ho change by more than 
1.5-2g . Thi s shift agre e s with the findings in iKessler et al.l 
(|2009|) and iBengocheal (|201ll ) and suggests, in our opinion, 
that the choice of the model used to fit the light curves can 
bias sensibly the cosmological results obtained. 



Effects of supernovae systematics and light curve fitters 

Uncertainties in the modelling of SNIa can affect the cos- 
mological parameters and associated errors extracted using 
a given dataset. To test the impact of systematics and light 
curve fitters on the results just presented, we re-analyse the 
ACDM cosmology with two different SNIa settings. First 
we use the same Union2 set, but neglecting the system- 
atic errors provided with the data; second we substitute 
this dataset with the Sloan Digital Sky Survey-I I Supernova 
Survey sample (SDSS SN. iKessler et al.ll2009l l. The sam- 



6.2 Curvature 

In this section we analyse the cosmological constraints ob- 
tained when the curvature is considered as a free parameter. 
The full list of parameters (mean plus the 68% confidence 
level) is shown in Table[5]for the combination of experiments 
listed in the header of the table. 

The blue shaded areas within dashed lines in panel a) 
of Figure [TT] and in Figure [TS] show the two-dimensional 
marginalised constraints in the Ode — fiivi and Ode — ^0 
planes, respectively, as obtained from CMB data alone. The 
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Table 5. Marginalised constraints on the cosmological parameters of the kACDM parameter space from the combination 
of probes listed in the header of the table. The quoted values are the same as in Table [3] 
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Figure 11. Two-dimensional marginalised constraints of the kACDM parameter space in the Ode — f^M plane. Panel b) zooms into 
panel a) in order to highlight the constraints obtained when combining CMB with the other datasets used in this work. The diagonal 
dotted line is for a flat Universe. Colour and line coding is the same as in Figure [9] 



plots show a very strong degeneracy between these param- 
eters. As stated in Section 16.11 the apparent size of the 
CMB acoustic peaks depends on their physical size and 
the angular diameter distance. Given that curvature den- 
sity is small and that it scales as the inverse square of the 
scale factor, it affects much more Da than the dynamics 
of the early universe. It is always possible to find combina- 
tions of fik, ijJm and ujh that keeps the angular size of the 
CMB peaks constant. Therefore, the constraints on the de- 
rived parameters are much weaker than the ones obtained 



in the previous section: ft-oE 
and Ho = 57^12 kms"^ Mpc~ 



fiiv 



The measured curvature 



0.451°:?° 



"k = (-4.9+t.t) 
the Icr level. 



X 10 ^ , is compatible with flatness at about 



Both panels in Figure [Til and Figure [121 show, with red 
shaded areas enclosed in solid lines, the same parameter 
planes as above when explored combining CMB and large 
scale structure information. Panel b) of Figure [TT] shows a 
zoom of the constraints of panel a) in the area around Ode = 
0.73 and Hm = 0.28 in order to show with more details the 
constraints obtained when more information is added to the 
CMB data. The galaxy power spectrum is very sensitive to 
the matter density. Additionally the BAOs allow to measure 
the angular diameter distance to the mean redshift of the 
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Figure 12. Two-dimensional marginalised constraints of the 
kACDM parameter space in the f^DE ~ Ho plane. Colour and 
line coding is the same as in Figure |9] 



LRG sample, z = 0.313. The degeneracies in the CMB alone 
are therefore strongly reduced when LSS measurements are 
included. The four variables discussed in the previous para- 
graphs become ^de = 0.731 ± 0.019, Qm = 0.267Io:o2i> 
Ho = 71.3 ± 2.5 and ^k = (1.6 ± 5.4) x 10"^: their errors 
are almost one order of magnitude smaller than in the CMB 
only case. With respect to the corresponding ACDM case, 
the uncertainties on the first three quantities increase up to 
50%. 

The inclusion of SNIa and Ho measurement, decreases 



+4.81 



xlO" 



the errors on curvature by about 10% (f2k ~ [3 
and on the other three parameters considered before by 
circa 20% (Ode = 0.7341°;°^^, Qm = 0.263to:ol5 and 
Ho = 71.9^20)- When comparing with the results from the 
previous section, the errors on dark energy and matter den- 
sity are unchanged, while increasing by less than 40% for 
the Hubble parameters. The two-dimensional marginalised 
constraints in the JIde — ^m and Qoe — Ho planes are shown 
in panel b) of Figure [TT] and in Figure [T^ with green shaded 
lines enclosed in dot-dashed lines. 



6.3 Beyond the cosmological constant 



In sections 16. II and 16. 21 we assume that dark energy is mod- 
elled as a cosmological constant. Despite its simplicity and 
its success in describing simultaneously many independent 
observations, the present day value of the density of dark 
energy does not have a solid physical explanation; this led, 
in the past decade, to the exploration of a large number 
of alternative model, most of which present a time depen- 
dent equation of statqj. Ideally, one would like to be able to 
constrain the full time, or redshift, dependence of woEiz) 
in order to restrict the range of possible models. Usually, 



° All theories of modified gravity can be also represented through 
an effective WD^i^)- 



parametric forms for the dark energy equation of state are 
assumed, which allow to measure time dependencies, but do 
not necessarily reproduce the correct u;de(z). With non- 
parametric approaches and principal component analysis 
it is possible to overcome some of th ese limitation s (e.g., 
Huterer fc Starkmanir2003l : ISerra et al]i2009 : Holsclaw et all 
20101 ). In this section we assume the simplest parametric 
form possible: kjde is a constant, independent of time, with 
a flat prior in the range [—2,0]. Deviations from the value 
Woe = — 1 would suggest that the cosmological constant 
is not a viable model of dark energy. If this is the case, 
the results shown in the previous two sections could be bi- 
ased as consequence of wrongly assuming that dark energy 
is described by A. In section 16.51 we analyse an alternative 
scenario in which iude is parametrized as a linear function 
of the scale factor. 

Table[H]lists the constraints on the cosmological param- 
eters as obtained for the different combinations of datasets 
analysed in this work. They are in perfect agreement with 
the ones presented with the previous two parameter spaces. 

Panel a) of Figure [13] shows the constraints in the 
Ode — wuE plane. The horizontal dotted line corresponds 
to iode ~ —1. If the CMB alone is considered a strong de- 
generacy between these two parameters is present, as shown 
by blue shaded areas enclosed by dashed lines. The reason 
for this degeneracy is analogous to the case analysed in the 
previous section when curvature is left as free parameter. In 
fact, while dark energy is subdominant at early times and 
does not influence the physical scale of the acoustic peaks, 
its equation of state parameter and density have a strong 
effect on the angular diameter distance to the last scatter- 
ing surface. Thus different combinations of lubm, i^b and 
wuE can result in the same apparent position of the acous- 
tic oscillations in the CMB temperature power spectrum. 
The constraints on the dark energy density and equation of 
state parameter are therefore very weak: ^de = 0.641q {g 
andiUDE = -0.74:tofo- 

With the inclusion of the LRGs data, the degeneracies 
in the CMB are broken, as shown in the panel a) of Figure 
[13] (solid line enclosing red areas). From the same figure it is 
clear that the impact of supernovae and Ho (dot-dashed line 
enclosing green areas) information is much larger in this sce- 
nario than in the previous two. In order to understand what 
causes such a difference, we show in panel b) of Figure [T51 also 
the two dimensional constraints obtained from the combina- 
tion of CMB+PLRG(fc)+SNIa (dashed lines enclosing blue 
areas). The horizontal line corresponds to woe = — 1. The 
plot shows clearly that supernovae have a large constrain- 
ing power for this parameter space. In fact the supernovae 
luminosity-distance relation traces the expansion history of 
the Universe and clearly identifies the transition betwee n 
the deceleration and acceleration phases (jRiess et al.ll2004l '). 
This transition depends on the densities of the cosmic com- 
ponents and is very sensitive to the dark energy equation 
of state. The inclusion of a precise measurement of Ho in- 
creases the precision on the parameters by a small factor. 
The constraints on the dark energy equation of state are 
i) wuE = -1.02 ± 0.13, from CMB4-PLRG(fc), ii) wde = 
-1.07±0.11, from CMB+Pi^nGik)+Ho, iii) wde = -1.009± 
0.069, CMB-hPLRG(fc)+SNIa and iv) woe = -1.025to:o65> 
from CMB+PLRG(fc)+-H'o+SNIa. Therefore the inclusion of 
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Table 6. Marginalised constraints on the cosmological parameters of the wCDM parameter space from the combination of 
probes listed in the header of the table. The quoted values are the same as in Table [3] 
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Figure 13. Panel a): two-dimensional marginalised constraints of the wCDM parameter space in the uide — f^DE plane. Colour and 
line coding is the same as in Figure |9] Panel b): same as the upper panel, but for CMB-l-PLj^G(fc), red areas within solid lines, 
CMB+PLjjQ(fc)+SNIa, blue areas within dashed lines, and CMB+PLptQ(fc)-|-iyo+SNIa, green areas within dot-dashed lines. Note that 
the lower panel is a zoom of the upper one in the region occupied by the CMB+PLp^Q(fc) contours. The dotted horizontal line shows 
w^DE = -1 



the large scale structure measurements decreases the error 
on WOE by about a factor 3 and SNIa further halve it. 

The main result of this section is that, when combin- 
ing CMB, LSS, SNIa information with the prior on Ho, the 
equation of state parameter of dark energy is constrained 
to be -1 with about 6.5% accuracy. This, although perfectly 
compatible with the cosmological constant models, does not 
exclude other dark energy scenarios. Excluding CMB, which 
has very little constraining power, the other three experi- 
ments give us access only to redshift z < 1. Furthermore 



many dark energy models can be tuned in order to mimic 
A at low redshifts and deviate from it at earlier epochs. In 
order to narrow the range of possible models, more precise 
measurements spanning a larger range of redshifts, together 
with the inclusion of tests of the growth of structures, will 
be necessary. 
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Table 7. Marginalised constraints on the cosmological param- 
eters of the wCDM parameter space from the combination 
of CMB+PLf^Q(fc)+i?o+SNIa when systematic errors are not 
considered and when the mlcs2k2 SNla light curve fitter is 
used. 
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Figure 14. Two-dimensional marginalised constraints of the 
kwCDM parameter space in the iude ~ f^k plane. Colour and 
line coding is the same as in Figure [9] The horizontal and verti- 
cal dotted lines show the values of these parameters for the flat 
ACDM case. 



Effects of supernovae systematics and light curve fitters 



Similarly to what has been done in section [6711 we test the 
impact of SNIa systematic effects and different light curve 
fitters on the measured cosmological parameters. The con- 
straints for these two cases are listed in Table [T] As be- 
fore, we find that neglecting systematics does not change 
the mean values of the parameters, but generally reduces 
the associated errors. In particular the constraints on the 
dark energy equation of state are reduced by almost 30%, 
wde = —1.007 ± 0.046. The use of the data from the 
SDSS SN project, with the light curves fitted with mlcs2k2, 
changes some of the parameters, as for instance Ho, u;de and 
Q,M , by more than 2-cr with respect to the value in our stan- 
dard case. In particular we measure wde to be — 0.875lo!o55- 
This highlights, better than for the ACDM case, the im- 
portance of SNIa modelling in improving cosmological con- 
straints from future generation experiments. 



6.4 Curvature and dark energy equation of state 
as free parameters 

We now analyse the accuracy that can be achieved when 
both wuE and Slk are considered as free parameters. The 
full list of constraints on the cosmological parameters for 
the kwCDM cosmology is summarised in Table [S] 

The introduction of an extra degree of freedom with 
respect to the previous two sections affects the CMB degen- 
eracy already discussed. The two-dimensional marginalised 
constraints in the wde — i^k and uide — S^de planes from 
CMB data only are shown in Figures [14] and [15] with the 
blue shaded areas within dashed lines. By comparing Figure 
1151 with panel a) of Figure [T3] it becomes clear that the al- 
lowed area in the parameter space increases because of the 
larger degeneracy. This causes the errors on fl^ and uide to 




DE 



Figure 15. Two-dimensional marginalised constraints of the 
kwCDM parameter space in the «>de — f^DE plane. Colour and 
line coding is the same as in Figure |9] The vertical dotted line 
shows «)DE = — 1 



increase by 50-60% with respect to the corresponding cases 
in sections 16.21 and 16.31 

The red contours within solid lines in Figures [14] and 
1151 show the constraints in the wde — ^k and woe — f^DE 
planes from the combination of CMB with the LRG data. 
The inclusion of the power spectrum reduces the region al- 
lowed by CMB alone to a one-dimensional degeneracy, that 
can be broken using the information on the amplitude of 
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Table 8. Marginalised constraints on the cosmological parameters of the kwCDM parameter space from the combination 
of probes listed in the header of the table. The quoted values arc the same as in Table [3] 
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PhRaik). Therefore a degeneracy arises between the bias 
and the shape of the power spectrum when both the cur- 
vature and the dark energy equation of state are treated 
as free parameters. When considering the range of scales 
0.02 /iMpc"^ 5? fc «; 0.15 /iMpc"^ larger values of woe, Om 
and f2k and lower values of as can be compensated by un- 
physically large biases 6 and larger values of Amc in order 
to obtain a comparable x^- H we used a finite flat prior on 
b we would decrease the large f^k and «jde tail and reduce 
the degeneracy. Alternatively, the detection of the turnover 
in the power spectrum, not possible nowadays because of 
the still too small volumes probed by galaxy surveys, might 
help constraining better its overall shape and break the de- 
generacy just described. 

The addition of Ho and supernovae measurements 
breaks the bias-shape degeneracy, returning parameters per- 
fectly consistent with the ACDM cosmology. As for the 
previous section, the biggest change in precision is due to 
the SN, which improves the accuracy of the three param- 
eters previously discussed, by a factor 2-3. The final con- 
straints, when all four probes are used, are shown in Fig- 
ures [T?] and [TJ] by the green shaded areas enclosed within 
dot-dashed lines. For the quantities shown in the plot we 
obtain fik = (4.5 ± 6.5) x 10""^ (36% increase in the errors 



with respect to kACDM case), «;de = — 0.981_ 



(29% 



increase in the errors with respect to the wCDM case) and 
J^DE = O.733I0017 (8% increase with respect to both). 



Time varying dark energy equation of state 
parameter 



6.5 



In section 16.31 we encode possible deviations from the cos- 
mological constant model in a constant effective dark energy 
equation of state. In this section we include explicitly a time 
dependency on woe through the simple parameterization of 
equation ^. We treat the two parameters of this model, wq 
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Figure 16. Two-dimensional marginalised constraints of the 
waCDM parameter space in the wg — Wa, plane. Colour and line 
coding is the same as in Figurc[9] The vertical and horizontal dot- 
ted lines show the values of these parameters for the fiat ACDM 
case (too = — 1 and tOa = 0). The dotted diagonal line shows the 
equation lODE(ap) = —0.97 = wo + (l — ap)ws., where ap = 0.65 is 
the pivot scale factor from the combination of CMB and PLRG(k) 
(see text). 



and Wa as free and we consider the curvature to be fixed 
0. The constraints on the cosmological parameters for the 
cosmology analysed in this section from the different com- 
binations of probes are listed in Table |9l 

Figure [T6l shows the two-dimensional marginalised con- 
straints in the wo — Wa plane from CMB alone (blue shaded 
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Table 9. Marginalised constraints on the cosmological parameters of the waCDM parameter space from the combination 
of probes listed in the header of the table. The quoted values are the same as in Table [S] 
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areas within dashed lines), CMB plus LRG power spec- 
trum (red shaded areas within solid lines) and the com- 
bination of CMB, PLRG(fc), ^0 and SNIa (green shaded 
areas within dot-dashed lines). The vertical and horizon- 
tal dotted lines at u)o = — 1 and Wa, = Q show the val- 
ues of the two parameters for the ACDM case. In the three 
cases a degeneracy is visible, which is reduced as more in- 
dependent data are included. From CMB alone we obtain 
Wo = -0.7llJJ;4g and Wa = -Q.2,2t\°o\; the inclusion of 
the LRGs power spectrum information does not change sub- 
stantially the one-dimensional constraints (tuo = — 0.75lo Is 
and Ws, = — 0.63lo 9g) but reduces by almost a factor 2 the 
figure of merit, i.e. t he area of the 95% confidence level 
(jAlbrecht et al.l 120061 ). As shown in the previous two sec- 
tions, supernovae have an important role in constraining 
dark energy properties and the inclusion of their infomation 
reduces the errors in the two parameters by about a factor 
3 and 2, respectively. The one-dimensional constraints from 
the combination of all the datasets are wo = —1.00 ± 0.14 
and Wb, = —0.13 ± 0.53. The attempt to constrain also the 
time evolution of the dark energy equation of state results 
in a degradation of its present value. 

From these results it is possible to reconstruct the time 
dependence of the dark energy equation of state parameter. 
The thick dashed, solid and dot-dashed lines in Figure 1171 
show the value of tode(2) from CMB, CMB plus Phucik) 
and the combination of the four probes, respectively. The 
corresponding 1-cr errors, which vary with redshift, are in- 
dicated by the shaded areas within the thin lines. They are 
computed according to lAIbrecht et al.l (|2006l ): 



((5tUDE(a)) = {{Swq + (1 - a)5wa) )• 



(16) 



In the redshift range shown in the plot, wue is al- 
ways compatible with a cosmological constant at the la 
level. Furthermore th e errors show a minimum at a red- 
shift called "pivot" (JHuterer fc Turneij 1200 ll : IHu fc JainI 
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Figure 17. uide(^) as function of the redshift z as obtained from 
the constraints shown in Figure [16] The thick dashed blue, solid 
red and green dot-dashed lines are the mean uideC^) for CMB, 
CMB-I-PlrgCc) and CMB-|-PLRG(fc)+fl'o+SNIa. The blue, red 
and green dashed lines enclosed within the outermost dashed, 
solid and dot-dashed lines are the 68% confidence level for the 
same three combinations as computed from equation l|16| l. 



l2004l : lAlbrecht eral] l2006l l. For the three cases shown 
in Figure 1171 we obtain a pivot redshift of Zp — 0.4 
(CMB), Zp = 0.54 (CMB-hPLRG(fc)) and Zp = 0.3 
(CMB-|-PLRG(fc)+^o+SNIa) and an equation of state pa- 
rameter woe{zp) = -0.80 ± 0.37, woEizp) = -0.97 ± 0.29 
and w-de{zp) — — 1.03 ± 0.069. It is interesting to note that 
in the last case the precision is comparable to the one pre- 
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Figure 18. Two-dimensional marginalised eonstraints of the 
waCDM parameter space in the uiq — f^DE plane. Colour and 
line coding is the same as in Figure |9] The vertical line is for 
wo = -1. 



sented in Section [6.31 The diagonal dotted line in Figure [161 
shows that the degeneracy between wo and w^ is very close 
to Wo + (1 — ap)ws. — —0.97, where Sp — 0.65 is the pivot 
scale factor from the CMB-l-PLRG(fc) constraints. 

Figure [18] shows the two-dimensional marginalised con- 
straints in the wg — Ode plane. The colour and line code 
is the same as in Figure 1161 This figure illustrates that the 
addition of large scale structure information to the CMB 
data halves the errors on the dark energy density from 
Qde = 0.66 ± 0.10 to Qde = 0.702tggg. The inclusion 
of SNIa and Ho measurements decreases the errors further 
by 70% to f^DE = 0.734 ± 0.016. 



7 COMPARISON WITH PREVIOUS STUDIES 

In this section we compare our results with recent work fo- 
cused on the analysis of the large scale structure of the Uni- 
verse. In Section lTTll we p erform a quantitat ive comparison of 
our results with those of lReid et al.l (|20ld . thereafter RIO), 
whose measurement and model are publicly availablqj. Sec- 
tion l7.2l instead, lists a number of similar analyses, for which 
this is not the case, and we can only discuss the diflerences 
qualitatively. 



7.1 iReid et al] (|2010l ') 

The LRG distribution from the SDSS DR7 has already been 
used to extract cosmological parameters by RIO. The anal- 
ysis in that work differs from the one we perform mostly in 
five important details; i) all the LRGs from the Northern 
Galactic Cap and the 3 southern stripes were used (110576 



galaxies in 7931 deg 



thanks to the count-in-cylinders 



'^ |http : //lambda . gsf c ■ nasa ■ gov/toolbox/lrgdr/ 1 
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technique (CiC. lReid fc SpergejIJOOgl : iReid et al.ll2009l ), the 
authors extracted a halo catalogue from the LRG distri- 
bution and then computed the halo power spectrum using 
the PVP estimator, iii) they computed the covariance ma- 
trix from 10000 lognormal catalogues (LN, IColes fc Joned 



1991 ). iv) they used a model based on halofit ( Smith et al.l 
20031 ) ■ that required additional calibration against numeri- 
cal simulations and v) they analysed the power spectrum 
in the range 0.02/iMpc~^ < A; < 0.2/iMpc'\ The ad- 
vantage of using the reconstructed underlying density field, 
instead of the galaxies, is that the intra-halo peculiar mo- 
tions, which cause the so called fingers-of-gods, are erased, 
leaving weaker small scale redshift-space distortions. They 
combine their halo power spectrum with the 5th year data 
from the WMAP satellite fWMAP5. iKomatsu etal] 120091 ') 
and the Union Supernovae dataset (JKowalski et al.l 120081 '). 
The comparison between the results in table 3 of RIO and 
the corresponding ones in this work (third column in Ta- 
bles |3l O |6] and third and fifth columns in Table [8]) shows 
that, with the exception of the ACDM case and despite the 
smaller k modes used by us, the errors that we obtain here 
are slightly smaller and that there are significant offsets in 
the preferred values of some parameters. For example, the 
constraints on wue in the wCDM parameter space, when 
combining CMB and P(k), are woe = —1.02 ± 0.13 in our 
analysis and woe = —0.79 ± 0.15 in RIO. 

As the full measurement from RIO and a COSMOMC 
module which includes the model is publicly available, we 
can understand the origin of the differences between the 
two studies. We concentrate here on the wCDM parame- 
ter space. Table [TOl compares the results shown in the third 
column of Table [6] (second column) with the constraints that 
we obtain applying i) RIO model to our measurement (third 
column), ii) our model to RIO measurement (fourth column) 
and iii) RIO model to RIO measurement (last two columns). 
To do it, we combine the large scale structure information 
with the CMB measurements presented in Section 14.11 we 
consider the scales 0.02/iMpc~^ ^ fci ^ 0.15/iMpc~^ and 
we compute the model for kj ^ 0.2 feMpc"^ before the con- 
volution with the window function. In the last column we use 
0.02 /iMpc"^ < fc, ^ 0.20 /iMpc"^ and kj sj 0.50/iMpc"\ 
as originally done in RIO. The similarity between the con- 
straints obtained in this last case and those of RIO shows 
that the effect of using 'WMAP7 plus small angular scale 
CMB data instead of WMAP5 is very small. If we con- 
sider the parameters constrained mainly by the CMB (ojb, 
i^DM, 6, T, Ab, 2ze and Us) we find that they are, as ex- 
pected, mostly independent of the power spectrum measure- 
ment and model used in the analysis. Excluding the last 
column, for the other parameters we find i) that, given a 
measurement, the errors that are obtained using the model 
of equation Q are generally smaller than the ones from 
RlO's model and ii) that, given a model, the parameters ob- 
tained from our measurement agree better with the bulk of 
the literature than the ones recovered from the RIO power 
spectrum. In particular, these data prefer larger values of 
wuE and of Qm, which lead to smaller values of ug and Hq. 
Using PLRG(fc) presented in Section lSTTI we measure the dark 
energy equation of state parameter to be 'W„„_ , „„„+o.i29 

when using our model and — 1.056lo'i9i when using RIO 
model. This is more than a 50% increase in the uncertainty. 
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Table 10. Comparison between this work and RIO: marginalised constraints on the cosmological parameters of the 
wCDM parameter space from the combination of CMB and P{k). The central four columns are for 0.02hMpc~^ ^ ki ^ 
0.15/iMpc-l and kj ^ 0.2/iMpc-^ while the last one for 0.02hMpc-l ^ fc; C 0.20 /iMpc"! and kj ^ O.SfeMpc^^ The 
quoted values are the same as in Table |3] 
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If we consider the RIO measurement, we obtain, for the 
two models, Wde = -0. 7981 jiJ: 127 and -0.847lo:i^g, with 
wbe > —1 at 1.5 and 0.9 a level, respectively. Including 
scales up to 0.2/iMpc~i, we measure wde = — 0.817lo:jg2: 
the error decreases, but remains larger than when using our 
model, and disfavours the cosmological constant scenario at 
1.2a. These results clearly point towards differences in both 
the modelling and the measurement of the power spectrum. 

While equation ((Sjl makes use of only three free pa- 
rameters and does not require any calibration, the model in 
RIO relies largely on numerical simulations in order to fix 
some of the parameters and to correct for the imprecision 
of HALOFIT, which at fc ^ 0.1 — 0.2/iMpc^i can describe 
a A CDM power spectrum with about a 5% accuracy (see 
e.g. iHeitmann et al.l 120101 ). Furthermore the calibration is 
performed only against simulations that reproduce the clus- 
tering of a flat ACDM universe, which could potentially in- 
troduce biases for more general parameter spaces, as the 
comparison between the second and third column in Ta- 
ble [TO] might suggest. The parameters shown in Tables \EIS\ 
are almost insensitive to the cosmological model assumed, 
whilst the ones listed in table 3 of RIO exhibit larger vari- 
ations. This might indicate that our model is better suited 
to analyse a large variety of cosmological parameter spaces. 

The tension between the constraints obtained when us- 
ing different measurements, instead, points in the direction 
of differences in the shape of the power spectrum or of the 
window matrix. Figure [19] shows the LRG power spectrum 
described in Section [3.11 (with blue points) and the halo one 
from RIO (red diamonds). The error-bars show the corre- 
sponding variances. Besides a factor of about four in the 
amplitude, due to the different samples and estimators used, 
which we include in the figure, the shapes of the two power 
spectra are almost identical, as we expect. We would also ex- 
pect the window functions to be similar, but Figure [2I]l shows 
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Figure 19. Comparison of the LRG power spectrum described 
in Section 13.11 (blue points) and the halo one from RlO (red di- 
amonds) with the corresponding errors. The amplitude of the 
power spectrum form RlO has been increased by a factor 4 in 
order to make the comparison easier. The best fit linear power 
spectrum of Figurc[2]convolvcd with the window matrix presented 
in Section 13.11 and from RIO are shown with a blue dashed and 
a red dot-dashed lines, respectively. In order to make the plot 
more readable, the amplitude of these two power spectra is not 
matched to the one of Plrg- 



that the rows of W{ki,kj) corresponding to the fc-bands of 
the measured power spectra are substantially different. This 
is intriguing, given that the bulk of our and RIO samples 
the same. This difference is even clearer when convolving a 
model P{k) with the window function. The dashed and dot- 
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Figure 20. Comparison of tlie rows of the window matrix corre- 
sponding to the fc-bands of the measured LRG power spectrum 
analysed in this work (solid lines) and from RIO (dashed lines). 
For clarity only one every fourth row is shown. 



dashed lines in Figure[T9]show the best fit linear power spec- 
trum from Figure [2] after the convolution with our and RIO 
W{ki,kj), respectively. In order to avoid clutter the ampli- 
tude of these two power spectra is not matched to PLRG(fc)- 
The differences in the window matrix are transformed into 
different shapes of the convolved Piin(fc), which implies that 
these measurements cannot be described by the same set of 
cosmological parameters. This difference in the window ma- 
trix is responsible for the main differences shown in Table 

m 



7.2 Other works 



IPercival et al.l (|201Cl ) analysed almost 900,000 galaxies from 
the combination of the full SDSS DR7 galaxy sample and 
the two-degree Fi eld Galaxy Redshift Survey (2dFGRS, 
ICoUess et al.l [2003J) and extracted the BAO feature from 
the power spectrum in 7 redshift bins. The main cosmo- 
logical results are shown in their table 5 and can be com- 
pared with the third column in Tables [S] [5] and [6] and with 
Table [8] Their constraints are compatible with ours and 
with the findings of RIO. Significantly, despite the much 
smaller sample that we use, the constraints that we obtain 
are co mparable or tighter than the ones of IPercival et al.l 
(|2010h . As an example, for the wCDM case we obtain from 
CMB-l-PLRG(fc) WDE = -1.02 ± 0.13, while they report the 
value Woe = —0.97 ± 0.17, which corresponds to an im- 
provement of about 25% in our results. This comparison 
suggests the importance of the use of the full information 
content in the galaxy po wer spectrum, as already noticed 
by different authors (e.g.. ISanchez. Baugh fc Anguld 120081 : 



IShoii. Jeong fc Komatsi3l2009l : iBlake et al.ll201ll ). 

The correlation functi on has also b een in tensively 
used with similar goals. I Sanchez et al.l ((20091) applied 
a model equivalent to the one presented in section 
S.llJCrocce fc Scoccimarrdl2008l : ISanchez. Baugh fc Anguld 



20081 ) to the correlation function of the LRG sample from the 



SDSS DR6, as measured by ICabre fc Gaztanagal (|2009al ). 
Combining this with Union SNIa sample and the CMB mea- 
surements from WMAP5, as well as with the position of the 
BAO peak along the line of sight in the two-dimensiona l 
correlation function from iGaztanaga. Cabre fc Huil ((20091), 
they extracted cosmological parameters for the same pa- 
ra meter spaces p r esente d in this work. The overall results 
of ISanchez et al.l ((20091 ) are consistent with the ones pre- 
sented in section [6| Interestingly our errors on the cosmo- 
logical parameters are systematically smaller when curva- 
ture is kept as free variable and larger for the wCDM and 
waCDM cases. This shows that the models and the mea- 
surements of the power spectrum and the correlation func- 
tion, although generally coherent, show some small differ- 
ences in the sensitivity to cosmological parameters for differ- 
en t parameter spaces. I n the kwCDM case, the constraints 
of [Sanchez et al.[ ((20091 ) show the same bias-shape degen- 
eracy that shifts the values of fik and wde upwards with 
respect to the fla t ACDM para d igm (compare Figure [T4[ 
with flgure 14 in (Sanchez et al.( ((2009( )). The inclusion of 
SNIa and of Ho (radial BAO) in our (their) analysis leads to 
the tightest constraints in both works. For flat and non flat 
ACDM, the differences are negligible, while for the other 
three cases, where SNIa measurements play an important 
role, the errors that we obtain are larger: this is due to 
the inclusion of th e supernovae s ystem atic errors, which 
where neglected in (Sanchez et al.l ((2009( ). in our analysis. 
When we do not use the systematics, we obtain tighter con- 
straints, shown by the middle column of Table [7[ In this 
case we measure the dark energy equation of state param- 
eter to be wde ~ —1.007 ± 0.046, ab out 11% tighte r than 
the corresponding value measured by (Sanchez et al.l ((2009( . 
WDE = -0.969 ±0.052). 

A more recent analysis of the correlation func- 
tion from the DR7 LRG samp le has been made by 
(Chuang. Wang fc Hemanthal ((2010h. using a sim plified ver- 
sion of the model presented bv (Reid et al.( ((20ld ). They com- 
bined this measurement with WMAP7 and Union2 data in 
order to extract cosmological parameters. Some of their con- 
straints are offset by l-2cr with respect to our findings and 
their measurements of the dark energy equation of state are 
compat ible with ours, once the SNIa systematics are ne- 
glected. (Carnero_e^_alJ ((201l( ) measured the angular correla- 
tion function and detected the BAO feature from the LRGs 
in the SDSS DR7 photometric catalogue, which consists of 
a sample of about 1.5 million galaxies. Fixing all the other 
parameters, with the exception of Hq, to the best fit from 
WMAP7, they measure tode to be -1.03 ± 0.16. 

The small scale clustering has Also been used to per- 
form cosmological analyses. From the small scale projected 
correlation function and mass-number r atio in cluster s, mod - 
elled in the halo model framework, (Tinker et al.l ((201 1( ) 
extracted cosmological and model parameters, obtaining, 
when combining with WMAP7, ^m = 0.290 ± 0.016 and 
(Tg — 0.826 ± 0.02. The difference between these results and 
the values of Table [3[ might be due to the smaller number of 
cosmological degrees of freedom and the higher number of 
model parameters. 

In a recent article (Blake et al.l ((201ll ) presented 
the first cosmological results from the WiggleZ survey 
((Drinkwater et al.l (20ld ). They used a sample of about 
130,000 emission line galaxies across 1000 deg^ in the red- 
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shift range 0.3 < z < 0.9. From the correlation function, the 
power sp ectrum, the BAO s and the band-filtered correlation 
function l|Xu et al.ll20ld l they extracted BAO parameters, 
like the effective distance Dy{z) of equation (|15|l and the 
"acoustic parameter" , defined by A{z) = Dv^O,m.Hq/cz 
l|Eisenstein et al.ll2005l V From those, they measured cosmo- 
logical parameters for the kACMD and the wCDM mod- 
els. When using LSS information only, the uncertainties are 
large but the results show a strong preference for the pres- 
ence of dark energy, with tode = — LG^o't- Combining the 
BAO parameters with WMAP7 distance priors they mea- 
sured the dark energy equation of state parameter to be 



wde 



-0.982: 



J, which became wde 



^1.026 ±0.081 



if also Union2 SN are used. Those results are in agreement 
with ours, although their uncertainties are about 20-30% 
larger than the ones that we show i n Table [6l The agre ement 
between our results and the ones in lBlake et al.l (|201in show 
that, despite the differences in the galaxies selected, the sur- 
vey volume and geometry and the procedure used to extract 
cosmological information, both analyses are robust enough 
for the precision achievable today. 



8 SUMMARY AND CONCLUSIONS 

We have computed the power spectrum of the distribution 
of about 90,000 luminous red galaxies, extracted from the 
spectroscopic part of the seventh data release of the Sloan 
Digital Sky survey. To compute the covariance matrix for 
the power spectrum we make use of the 160 LasDamas mock 
catalogues, which have the same angular and redshift dis- 
tribution as the observed LRGs. We describe the measured 
power spectrum with a model inspired by renormalised per- 
turbation theory and modified in order to describe biased 
objects in redshift space (section [STT|. This model has been 
successfully tested against numerical simulations (MIO) and 
mock catalogues (section l5.4p for k < 0.15 hMpc~^ at z=0- 
0.5. 

Combining the large scale structure information with 
measurements of the CMB temperature and polarisation 
power spectrum from the seven year data release of WMAP, 
ACBAR, BOOMERanC, CBI and QUAD and with the 
luminosity-distance relation from the Union2 supernovae 
sample and using a precise determination of the local Hub- 
ble parameter as a gaussian prior, we explore the constraints 
in five different cosmological parameter spaces described in 
section [S31 They are the fiat ACDM concordance model, a 
similar one where curvature is a free parameter (kACDM) 
and three models in which the dark energy equation of state 
has a parametric form: in two cases it is assumed to be con- 
stant, with and without the assumption of a fiat geometry 
(wCDM and kwCDM), and in the last case (waCDM) we 
model wde with the simple parametric formula of equation 

Overall, we obtain tight constraints on the cosmological 
parameters for all the five cases and we do not detect devi- 
ations from the fiat ACDM paradigm. The different combi- 
nations of the four experiments used in this analysis do not 
show any evidence of tensions between cosmological probes. 
We find that the curvature is null at l-cr level with errors 
of the order of 10~^ — 10""^. We measure the dark energy 
equation of state parameter to be consistent with a cosmo- 



logical constant with 13% uncertainty for CMB-f Plrg (fc) 
and 6.5% uncertainty for CMB-|-PLRG(fc)-l-//o-|-SNIa in the 
fiat wCDM case. If we discard the systematic errors in the 
SNIa, the precision increases to 4.6%. In the kwCDM, be- 
cause of the added degree of freedom, we have a degra- 
dation of the constraints to 8.4%, when combining all the 
four samples. The constraints obtained with the CMB and 
large scale structure together are shifted by 1.5-2a with re- 
spect to the best fit ACDM because of a bias-shape de- 
generacy in the power spectrum that allows very large, 
and unphysical, values of the bias when wue, ^m, ^k in- 
crease and as decreases. If we assume the parametric form 
of equation ^ for the dark energy equation of state, we 
obtain w-d-e{z = 0.54) = -0.97±0.29 (CMB-fPLRG(fc)) and 
wue{z — 0.3) = —1.03 ± 0.069 (all four experiments com- 
bined) . The latter is only slightly worse than the fiat wCDM 
result. 

In the near future new and larger galaxy redshift cat- 
alogues both spectroscopic, like BOSS and HETDEX, and 
photometric, like Pan-STARRS and DES, will become avail- 
able, together with the new measurem ents of the CMB 
anisotropies from the Planck satellite (jAde et al.l 1201 ll ). 
These datasets will enable us to improve the constraints 
presented in this work even further. Some of these experi- 
ments are explicitly designed in order to extract the max- 
imum amount of information regarding dark energy. This 
would allow to exclude many classes of dark energy mod- 
els or of modifications of general relativity. In order to use 
the full information from the power of these new large scale 
observations and to avoid introducing systematic effects, ac- 
curate models of the large scale distribution, scale dependent 
bias and redshift space distortions are necessary. In MIO we 
showed that the model used in this work is accurate enough 
to describe the power spectrum shape also for surveys with 
volumes larger than available nowadays. The explicit inclu- 
sion of bias and redshift space distortions in the model would 
allow however to use a larger range of scales than now possi- 
ble, helping to improve the quality of the cosmological con- 
strains even further. 
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APPENDIX A: BASIC EQUATIONS TO 
COMPUTE THE POWER SPECTRUM 

In this section we summarise the basic equations of the 
FKP and PVP estimators used to compute the power spec- 
trum and the window function from galaxy surveys. PVP is 
a generalisation of FKP and takes in account the change 
of the g alaxy bias with the luminosity of the g a laxies 



or trie g alaxy bias with the luminosity oi the g a laxies 
(see e.g.. lOavis fc Gelleil ll976l:jNorberg et all l200l( . I2OO2I : 



IZehavi et al.ll2002l : |Phleps et al.lf2006l ') 



The observed power spectrum Poik) is obtained from 
the squared average Fourier transform of the weighted den- 
sity field defined by 

1 



FKP: F(x) 



N 



w{x) [ng(x) -anr(x)]. 



(Ala) 



PVP: F(x) 



1 

iV 



dL 



w(x, L) 



.(x,L) 



r(x,L)], 



(Alb) 

where ng(x, L) and nr(x, L) are the number density of galax- 
ies and randoms of luminosity L at position x. The corre- 
sponding quantities of equation (|Ala|l can be obtained in- 
tegrating over the luminosity. w(x), iu(x, L) are weighting 
functions and 6(x, L), only PVP, is the bias relative to a 
specific galaxy population with luminosity L*. The normal- 
isation A'^ is defined by 



PVP:iV' = 



I d X / dLn(x, L)ui(x, L) 



(A2a) 



(A2b) 



where n(x, L) and ri(x) are, respectively, the mean expected 
number density, i.e. in absence of clustering, of galaxies of 
luminosity L at position x and its integral over L. Finally a 
is a constant introduced to match the two catalogues and is 
chosen requiring thai|f| {F{x)) = 0: 



FKP: a 



J d''a::ui(x)ng(x) 



PVP: a 



J d^x dL [to(x, i)/fe(x, L)] ng(x, L) 
J d^x dL [to(x, L)/b{x, L)] n,(x, L) ' 



(A3a) 



(A3b) 



The observed power spectrum can be then written for 
both estimators as: 

^o(k) = I |^P,(k')G^(k-k') = (|F(k)i^) - n„, (A4) 

where -Pt(k') is the "true" underlying power spectrum, the 
shot noise Psn is given by 



FKP: Rn = 



1 + a 



^— / d xn{'x.)w (x), 
J 



(A5a) 



In FKP, the authors use a definition of a such that A^ = 1. 
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PVP: Pen = 



1 + a 
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w2(x,L) 



25 



d xdLn(x, L) 



b\^,L) 



(A5b) 



and G (k) is the window function, which encodes informa- 
tion about the survey geometry. 

It can be shown that the window function is computed 
from the spherical averaged Fourier transform of the field 



FKP:G(x) = -n(x)«;(x), 



(A6a) 



PVP: G(x) = ^ I AL n(x, L)w(x, L), (A6b) 

as G'^{k) — (|G(fc)p) — Gsn- Gsn is the shot noise defined by 



FKP: Gsn 



1 

w 



d xn{x)w (x), 



PVP: Gs, 



1 

iV2 



d xdLn{x,L)w (x,L). 



(A7a) 
(A7b) 



In this work, we use the weighting function s, designed 
to min imise the variance, proposed in FKP and lCole et al.1 
(l2005fl : 



FKP:u;(x) = 



PVP:w;(x,L) 



1 + P(fc)n(x) 
Wib'^{x,L) 



1 + P{k) J dL 62 (x, i)n(x, L) ' 



(A8a) 



(A8b) 



The intrinsic weight of the objects, Wi, can con tain informa- 
tion about completeness and/or fibre collision (jZehavi et al.1 
[2002; Masicdi 20061'). P{k) is an estimate of the recovered 
power spectrum and it is usually substituted with a con- 
stant Pw, chosen in order to minimise the variance around 
the wave- number k for which P(k) ~ Pw 




p„=0, «;; = ! 

p„ =40000, «;; = ! 
p„ =40000, «;i=/c 
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APPENDIX B: IMPACT OF WEIGHTS ON 
POWER SPECTRA AND COSMOLOGICAL 
CONSTRAINTS 

In the following we test the impact of pw, Wi and estimator 
on the LRG and mock power spectra, on the errors and on 
the cosmological parameters. 



Figure Bl. Comparison between power spectra computed with 
Pw = 40000,0 (red dashed and blue sohd lines) and lo; = 1,/c 
(diamonds and up triangles) divided by a linear power spectrum 
without BAOs. The upper panel is for the FKP, the lower for 
PVP. The shaded area denotes the standard deviation computed 
from the mock catalogues using pw = 40000. 



Bl Testing the luminous red galaxies power 
spectrum 

In this appendix we test the impact of different choices of the 
estimator, pw and Wi on the LRG power spectrum and win- 
dow function. For both estimators, FKP and PVP, we test 
four values of Pw = 40000, 10000, 4000, 0. We also consider 
four different intrinsic weights: i) Wi = 1 (all the objects have 
equal weight), ii) Wi = c (areas with low completeness have 
less weight than areas with higher one), iii) Wi = fc (the loss 
of galaxies due to fibre collisions is compensated as described 
in Section [3.1(1 and iv) Wi = c x fc (both completeness and 
fibre collision corrections applied). 

Figure IBll shows the differences in the power spectra 
measured with the FKP (upper panel) and PVP (lower 
panel) estimators for the different choices of pw and Wi. 
For clarity, in the figure we show only the combination of 



Pw = 40000, (red dashed and blue solid lines respectively) 
and Wi — l,fc (diamonds and up triangles, respectively). 
The shaded area shows the standard deviation as measured 
from the mock catalogues for pw = 40000. All the power 
spectra have been divided by a non-wiggle one with the same 
cosmological parameters as the mock catalogues. Different 
choices of pw change marginally the shape of the power spec- 
trum. The results for pw = 10000, 4000 fall in between the 
two extreme cases shown in Figure iBll 

Instead, the correction for fibre collision has a siginfi- 
cant effect. When obtaining spectra of crowded fields, not all 
the objects of interest can be targeted with a limited number 
of pointings. Because of this loss of objects, the amplitude of 
the highest peak in the density fields decreases, while the low 
density regions are unaffected. This causes the amplitude of 
the fluctuations, and consequently of the power spectrum, 
to be lower. We have indeed measured a few percent scale 
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0.05 0.10 0.15 0.20 
k[h/Mpc] 

Figure B2. Ratio between power spectra computed the FKP and 
PVP estimators for pw = 40000, uii = fc. Tlie sliadod area shows 
the standard deviation computed from the mock catalogues for 
Pw = 40000. The dotted horizontal line is the mean of the ratio in 
the plotted interval. The other choices of Pf, and «); show similar 
behaviour. 

independent decrease in the amplitude of (|-F'(k)p) (i.e. the 
power spectrum before subtracting the shot noise) in the 
case when fibre collision correction is not applied. This dif- 
ference is visible at the large scale, low k, limit in Figure [Bl] 
where the power spectrum with the fibre collision (triangles) 
is always larger than the one without (diamonds). On the 
other side the shot noise f equations I A5|) depends on the ex- 
pected non-clustered number density (n) and the associated 
weights, which are not influenced (or influenced in a uniform 
way) by fibre collisions. This causes the shot noise to be the 
same in both cases, changing the shape of the final power 
spectrum Po(fc). This effect is clearly visible in the small 
scale, large k, limit in Figure IBf I where the fibre collision 
corrected power spectra become increasingly larger than the 
non corrected ones. We do not show the results when com- 
pleteness correction is included, since the measured power 
spectra overlap almost exactly the non corrected ones. This 
is because both the galaxy and the random catalogues are 
weighted in the same way, therefore the density fluctuation 
field -F'(k) is unchanged. 

Comparing the two panels of Figure [Bfl is it clear that 
the choice of the estimator has a strong impact on the re- 
covered power spectra: in the PVP case the amplitude is 
systematically lower than for FKP. Figure [B2l shows the ra- 
tio of the power spectrum computed with the latter estima- 
tor with respect to the one computed with the former for 
the case Pw = 40000, w^ = fc. As before, the shaded area 
corresponds to the standard deviation computed from the 
LasDamas catalogues. Although the amplitude is different, 
the relative bias between the two estimators is scale invari- 
ant in the range 0.02/iMpc~^ ^ fc sj 0.2/iMpc"^ This is 
expected from the fact that the galaxy sample that we use 
is almost volume limited and contains a very uniform pop- 
ulation of galaxies. 




10" 
k[h/Mpc] 



Figure B3. Window functions G^(fc) multiplied by the cube of 
the wave-number k as function of k for pw = 0, 40000 (solid and 
dashed lines, respectively) for ui; = 1. The upper panel is for FKP, 
the lower for PVP. 



Recently, iBalaguera-Antolmez et al.l (|20f0|) have shown 
that, in volume limited mock catalogues of the REFLEX ff 
cluster survey, the power spectrum computed with the PVP 
estimator has higher correlations than the FKP one already 
at fc > 0.f5/iMpc~^. Because of this and of the scale inde- 
pendent relative bias, we can safely use the power spectra 
as estimated with FKP in order to constrain cosmological 
parameters. 

The window function, describing only the radial and 
angular selection function of the survey, is not affected by 
fibre collision effect; on the other hand different choices of pw 
and the completeness correction change the weights of equa- 
tions (|A8|I . which influence the effective survey volume. As 
for the power spectrum, we do not measure differences when 
the completeness weighting is applied. Figure lB3l shows the 
"dimensionless" window function k^G^(k) computed with 
FKP (upper panel) and PVP (lower panel) forpw = 40000, 
(dashed and solid line respectively). Although the overall 
shape is similar, there are small differences in the oscilla- 
tions due to different effective volumes in the two cases. The 
results for pw = f 0000, 4000 fall in between the two extreme 
cases. The difference between the FKP and PVP window 
functions are negligible. 
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Figure B4. Comparison between mean mock power spectra com- 
puted for Pw = 40000, (red dasiied witii triangles and blue solid 
lines with diamonds respectively) divided by a power spectrum 
without oscillations. The shaded area denotes the standard devi- 
ation for Pw = 40000. 



When convolving these window functions with a hnear 
power spectrum as in equation ([l]), the differences in the 
range 0.02 /iMpc"^ < fc s^ 0.2/iMpc"^ are negligible. 



77 0.15 




1.0 



0.8 



0.6 



0.05 0.10 0.15 0.20 
k [h/Mpc] 



I 
I 



c 
o 



0.4^ 

o 
u 



0.2 



0.0 



B2 Testing the power spectra of the mock 
catalogues 

In order to test the impact of pw on the results from the 
LasDamas, we compute the power spectra of the 160 mocks, 
their mean, standard deviation and covariance matrix for 
Pw = 40000, 10000, 4000, 0. Figure [H shows the mean power 
spectra computed for the two extreme cases: the only dif- 
ference is a small change in amplitude, which confirms the 
findings of the previous appendix, i.e. that the impact of dif- 
ferent Pw on the computed power spectrum is negligible. The 
shaded area denotes the standard deviation for pw = 40000. 
The choice of p„ has, however, a large effect on the 
errors. Figures IB5I and IB6I show, respectively, the correla- 
tion matrix for pw = 40000 (upper panel) and pw = 
(lower panel) and their difference. At fc ~ O.lfeMpc"^, for 
Pw ~ 40000 we measure that the correlation is systemati- 
cally, although not significantly, lower than for p„ = 0. This 
is expected since, when using FKP, P{k ~ 0.1) ~ 40000. 
At larger wave-number the power spectrum amplitude is 
smaller than at fc ~ 0.1/iMpc~^ and the variance and the 
correlation are smaller for small values of pw- This justifies 
our choice of using pw = 40000 in our analysis, since the am- 
plitude of the power spectrum is of this order in the range 
of scales we are interested in. If the analysis were centred 
on the small scale power spectrum or correlation function, 
a smaller value of pw would be preferable. 

B3 Impact on the cosmological parameters 

Here we test how the differences in the shape of the LRG 
power spectrum and the covariance matrix influence the cos- 



Figure B5. Correlation matrices for pw = 40000, (upper and 
lower panels respectively). 
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Figure B6. Differences between the pw = 40000 and the pw = 
correlation matrices. 



mological parameters measured assuming the wCDM cos- 
mology. We combine the LSS information with the WMAP7 
data. 

Contrary to our expectations based on the results shown 
in Appendix IBll the cosmological parameters are much more 
sensible to pw than to Wi. For a fixed pw the covariance ma- 
trix is the same and only the shape of the power spectrum 
changes. As the shape of the power spectrum is not affected 
by the completeness correction, also the cosmological con- 
straints are insensitive to it. The fibre collision correction 
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instead changes the bias, over which we marginalise ana- 
lytically, without affecting the shot noise amplitude. Dif- 
ferences in the relative amplitude of the latter term can 
be absorbed at least partially by the mode coupling am- 
plitude ^MC, which is systematically, although not signifi- 
cantly, larger when the loss of galaxies due to fibre collisions 
is corrected for. Because of this, cosmological constraints 
remain almost unchanged for different wi . Changes of pw in- 
stead infiuence the covariance but only marginally the power 
spectrum. The differences in the former influence the cosmo- 
logical parameters, which differ by 0.5 — la, for pw = 40000 
and 0. For example I obtain the dark energy equation of 
state parameter to be iude = —1.02 ± 0.13 for pw = 40000 
and k;de = —1.10 ± 0.14 for p„ = 0. 

Given the precision that is possible to achieve with the 
data used in this article, the differences just highlighted are 
not distinguishable from the uncertainties in the parame- 
ters. But in future, given the big improvements expected, 
these effects might become important and will need further 
analysis. 

This paper has been typeset from a TJnjX/ IMJ5X file prepared 
by the author. 
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